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The  lowest  absorption  at  10. 6  p  found  with  GaAs  grown  from  the  melt  has  temperature 
dependence  expected  for  multi  phonon  absorption.  GaAs  made  by  chemical  transport 
showed  a  high  10. 6y  absorption  lue  to  inclusions.  The  previously  reported  value 
for  the  dielectric  constant  of  Gn\s  has  to  be  revised  to  E  =  13.6.  Diffuse  streaks 
in  the  diffraction  pattern  of  GaAs  was  attributed  to  diffuse  thermal  scattering 
of  electrons  at  acoustic  phonons.  Slip  in  GaAs  occurs  on  {111}  <110>.  Hardness 
of  GaAs  is  increased  by  tetragonal  defects. 


A  method  to  determine  oxygen  in  GaAs  with  a  sensitivity  of  5  x  1016cm  3  or  0.2  to 
1  x  a  monatomic  layer  of  0  was  developed. 


In  GaAs  ioped  with  silicon,  0  present  is  a  concentration  <  2  *  that  of  Si,  cannot 
be  determined  due  to  the  stability  of  Si02» 


The  theory  for  thermal  lensing  is  presented  in  modified  form.  A  method  for  the 
determination  of  mechanism  of  multiphonon  absorption  is  proposed. 
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ABSTRACT 


The  theory  of  thermal  lenalng  la  reviewed  and  presented  from  a 
different  viewpoint.  A  method  to  determine  the  contrlbntions  of  various 
processes  to  multiphonon  absorption  is  proposed. 

For  GaAs ,  the  lowest  10. 6p  absorptions  found  have  a  weak  temperature 
dependence  corresponding  to  theoretical  predictions  for  multiphonon 
absorption.  However,  other  mechanisms  might  give  the  same  dependence. 

Layers  of  GaAs  up  to  2mm  thick  have  been  made  by  chemical  vapor  transport. 

The  layers  had  inclusions,  however,  and  the  high  10. 6p  absorption 

(4  x  to  2  cm-i  Qr  larger)  probably  result  from  these>  Uslng  a  ngwly  developgd 

method  for  oxygen  detection  with  a  sensitivity  limit  of  5  x  1016cm-3 
or  0.2  to  1  x  a  monolayer  of  oxygen  show  ~  1018cm~3  0  in  GaAs-Te.  In 
GaAs-Si ,  an  oxygen  content  below  the  Si  content  -annot  be  determined. 

Repeated  dielectric  constant  measurements  on  GaAs  give  e  =  13.6.  Diffuse 
streaks  in  the  diffraction  pattern  of  GaAs  are  attributed  to  thermal  diffuse 
scattering  of  electrons  at  acoustic  phonos .  Slip  in  GaAs  occurs  on  {111} 
<110>.  The  hardness  of  GaAs  is  increased  by  tetragonal  defects.  Hot  pressing 
of  CdTe  in  presence  of  Cd(t)  still  gives  problems. 

KBr  vith  a  10. 6y  absorption  of  1.1  x  lO^cnT1  has  been  grown. 

A  TEA  laser  for  damage  studies  is  being  assembled. 
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1. 


INTRODUCTION 


j’hc  formal  of  this  report  follows 
report  in  which  projects  are  identified 
work  statement. 


closely  that  of  the  first  quarterly 
by  codes  keyed  to  the  contract 


Fhe  various  categories  are  briefly; 

a)  Crystal  growth 

b)  Polycrystalline  window  fabrication 

c)  Mechanical  properties  of  window  materials 

d)  Window  material  defect  characterization 

e)  Theory  of  residual  IR  optical  absorption 

f)  Absorption  measurement  techniques 

g)  Theoretical  evaluation  of  optical  performances  of 
windows 


3 


a- 1  Concentrations  of  Oxygen  In  Selected  GaAs  Samples 
P.  S.  Vijavkumar  and  J.  M.  Whelan 


The  previous  report  described  the  determintaticn  of  the  oxygen 
concentration  -  fugacity  relationship  in  Ga  solutions  for  oxygen  atomic 
fractions  below  1CT7.  This  information  and  the  measured  changes  in  oxygen 
fugacities  before  and  after  dissolving  known  amounts  of  GaAs  in  Ga 
solutions  have  been  combined  to  estimate  the  oxygen  contents  of  selected 
GaAs  samples.  The  apparatus  has  also  been  previously  described.  Fugacities  are 
determined  from  the  emf's  of  reverisble  (or  nearly  so)  cells  using  calcia 
stabilized  zirconia  as  the  oxygen  conducting  electrolyte.  Since  the  oxygen 

fugacities  are  small,  they  may  be  described  in  terms  of  the  equivalent  moleculr.r 
oxygen  pressures,  P 


0 


2  in  equilibrium  with  the  Ga  solution  at  the  zirconia 


surface  and  the  reference  oxygen  pressure,  PQ  in  equilibrating  HO  vapor  at 
-  •  -  2  ^ 


a  partial  pressure  corresponding  to  its  vapor  pressure  at  0°C  with  H2  at 
one  atomosphere,  i.e.  6  *  10~3  atm.  The  relationship  between  the  emf  (assumed 
to  be  reversible)  and  these  pressures  is 


1/2 


entf  =  |j~  *n| 


0„R 


Alternately  the  oxygen  fugacity  in  the  Ga  may  be  related  to  the  equivalent 
partial  pressure  of  H20  in  1  atom  of  H  : 

1/2 


ph2o  =  6  x  10 


A  GaAs  wafer  which  had  a  portion  of  its  surface  implanted  with  high  energy 


oxygen  ions  was  selected  as  a  standard  for  checking  the  accuracy  of  the  method. 


Figure  A  shows  the  variation  with  time  of  (P^/Pq  R)  1/2  before  and  after 


immersing  the  implanted  sample  in  the  Ga  solution  at  t  =  0.  Transit  effects 
associated  with  oxygen  gradients  in  the  Ga  solutions  upon  dissolution  of  the 
sample  were  minimized  by  extrapolation  of  the  long  term  drift  back  to  the  time 
the  sample  was  immersed.  Using  this  extrapolation  figure,  the  amount  of  oxygen 
in  the  melt  was  6  *  10 15  atoms.  This  is  to  be  compared  with  3  x  ic’5  atoms 
added,  estimated  from  the  amount  of  implanted  GaAs  dissolved  by  the  Ga 
solution.  Unfortunately,  the  latter  estimate  was  complicated  by  the  dissolution 
of  the  implanted  region  occuring  only  at  its  edges.  The  majority  of  the 
implanted  surface  was  not  wetted  by  the  Ga.  Several  explanations  have  been 


■*v»'&ytyi!?*slK%  "y  ip  --  '  ■*  ■  *' -  ! 


!c'  1?  U'W^-'V'  1 1**lV  r i" 


advanced,  the  ones  considered  mor L  likely  are: 

1)  the  implanted  surface  was  coveted  by  1  fllm  of  partlally  cracteJ/ 
polymerized  hydrocarbons  during  the  implantation 

2)  the  energy  of  the  implanted  ions  caused  the  formation  of  an  oxide 

such  as  gallium  arsenite  or  arsenate  which  is  not  wetted  nor 
dissolved  by  gallium, 

^  “f  °f  the  dlfflCUi“-’  agreement  to  within  a  factor  of  two  is 
ons  ere  satisfactory.  We  shall  continue  attempts  to  secure  reliable 
primary  oxygen  standards. 

Following  the  above  procedure,  a  heavily  Te-d„ped  GaAs  crystal  was 

InTn  1  v  r-  7  —  "  «■  —  -  be  between  0.3 

and  1  10  atoms  cm  3 ,  a  figure  grossly  in  excess  of  H  0 

es t ima ted  detection  limit  of  our  method  of  3  .  10-  oxygen  atoms  cm-. 

ternatively,  the  detection  limit  may  be  stated  as  a  fraction  of  0  2  -  1 

mono  layers  on  the  surface.  Dif ferentation  between  surface  and  bulb  oxygen 

an  be  made  in  principle  by  limiting  dissolution  to  a  shallow  surface  Zyer 
by  using  a  Ga  melt  ...  _  iayer 


v  w  - s-'-'  ^  biidiiow  surface  1; 

UT  3  Ga  "elt  a2m°st  saturated  with  GaAs.  figure  B  shows  the  time 

GsaTiTV'  „  °2/ V  bef0re  and  »fter  dlPPl"g  seiected  samples  of 

“  "e  “°UntS  °'  """"  -  «  the  detection  limits, 

Ifie  S 3TT1D I  P  nnrroonor,J4 _ _ 


Lne  aetectli 

The  sample  corresponding  to  curves  A  and  A'  were  those  for  GaAs 

Zth'i  prepared  by  Dr.  Jerry  „f  ^  ^  ^  ^ 

growth  in  a  auarr?  Kmot-  ^  y 


„T h  was  made  by 

*  3  qUartZ  b°at  "”d“  a  °V>  p«»«.  greater  than  or  eq„al  to 

'  pressure  of  Ga.O  in  equilibrium  with  0,0  at  900’C1.  This  sample 

as  previously  k„„„„  to  have  had  an  unusuaUy  low  trap  concentration2 
turve  B  was  nKtatnoa  .t  „  . 


f  - — w  tidp  concentration 

crys t  7  US1"8  thS  GaAS  Sa"Ple'  WA-10°°'  «»**.<  OaAs 

at  USC  7  8r°“n  f  qUart2  “Ith”Ut  dellbaratP  <*»W™  doping  by  W.  Allred 
SC  These  samples  were  amongst  the  ones  showing  the  lowest  optical 
absorption  in  the  10.6pm  region. 

Several  conclusions  were  drawn  from  the  above  results- 
(U  Gross  amounts  of  oxide,  i.e.  a  few  monolayers  as  a  possible  oxide 
residue  on  the  surface  of  GaAs  cleaned  by  acid  etches  and  water 
rinaes  are  removed  by  heatine  -In  H  -mm*.  *  j 
about  30  minutes  at  600°C  Thi  2  a  n  ng  a,  ppm  of  H  O  for 

the  .  r  atl"8  SChedule  18  that  to  equilibrate 

samples  before  partially  dissolving  them  in  Ga. 

“aS  det6Cte''  ln  3  8a">’le  «"»■  •»>«  circumstances  in 
ha  deliberate  attempt  was  made  to  provide  an  oxygen  overpressure. 

and  »7~::rr that  oxy8en  °ight  be  ~  ^ 

mpies  and  that  it  compensated  the  Si  in  these.  Under  these 
c  rcumstances  the  dissolution  of  Si  with  oxygen  in  gallium  might  lead  to  the 


■v»v.A-  ^ vb-s ct  -  a., . a.: -  * 


r 


IflShWl^  i,i  .  w  a 


formation  of  Si02  and  not  influence  the  emf  of  the  cell.  The  reaction  of 
Si  and  0  dissolved  in  Ga ,  based  on  the  typical  fugacities  of  oxygen  used 
and  the  assumption  of  ideality  for  Si,  should  proceed  to  completion.  As 
such,  oxygen  in  excess  of  twice  that  of  Si  would  only  be  detected.  The 
completeness  of  the  Si-0  reaction  also  suggested  that  the  method  could  be 
used  to  estimate  silicon  uncompensated  by  oxygen  in  GaAs . 

To  test  the  reason  of  Si  with  oxygen  in  Ga,  a  sample  of  Si-doped 
GaAs  was  dissolved  in  Ga  containing  a  known  amount  of  oxygen  as  determined 

by  the  cell  emf.  Evidence  for  the  reaction  of  Si  with  oxygen  is  shown 
in  Fig.  C. 

The  rapid  decrease  in  E  at  t  =  0,  the  conversion  time,  indicated  a 
decrease  of  oxygen  which  is  attributed  to  the  reaction  of  Si  and  oxygen 
and  transit  effects  associated  with  incomplete  mixing.  The  approach  to 
a  more  uniform  melt  is  associated  with  a  gradual  rise  followed  by  the  long 
term  drift.  Extrapolation  of  the  long  term  drift  assuming  the  Si-0  reaction 
product  to  be  Si02(s)  leads  to  an  estimate  of  the  Si  concentration  as  being 
4  x  ®cm-3 .  This  is  probably  in  fair  agreement  with  the  total  silicon 
concentration  of  the  sample.  Here  again  we  are  forced  to  draw  only  semi- 
quantitative  conclusions  because  of  uncertainties  in  the  standards. 

In  light  of  the  previous  results  it  is  not  surprising  that  we  were 
unable  to  detect  the  presence  of  a  suspected  oxide  caused  by  heating  GaAs 

in  air  to  'v  350°C  for  30  minutes,  a  procedure  known  to  decrease  the  optical 
transmission  at  10.6ym. 

References 

1.  Dr.  Jerry  Woodall,  private  communication 

2.  Dr.  Gordon  Kino,  private  communication 


Expitaxial  Growth  of  Semi-Insulating  GaAs  from  Ga  Solutions 
— *  -  T*  Hoelke>  P.S.  Vij  avkumar  and  J.  M.  Whelan 

A  horizontal  epitaxial  slider  unit  that  can  utilize  oxygen  pumping 

in  the  course  of  epitaxial  growth  has  been  designed  and  built.  This  and 

another  unit  will  be  available  shortly  for  the  growth  of  thick  films  which 

can  be  doped  with  Cr,  Fe  or  Ni  to  make  them  semiinsulating .  The  purpose  of 

growing  such  films  is  to  prepare  the  semi-insulating  crystals  in  a  manner  most 

likely  to  insure  selective  doping  and  its  effects  upon  optical  absorption  at 

10 . 6ym. 
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Figure 


-  2000  0  2000 
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Figure 


P.  J.  Sclicbta.  R.  E .  Chane 


a. 3a  Summary  of  Results  and  Plans:  Calorimetric  measurements  of  high- 
purity  melt-grown  KBr  cyrstals  indicate  a  10.6  micron  bulk  absorption 
coefficient  of  0.0011cm  which  is  still  much  higher  than  theoretically  pre¬ 
dicted.  In  order  to  more  surely  eliminate  contamination  by  bromate  ions, 
attempts  were  made  to  grow  crystals  from  molten  KBr  purified  by  reaction  with 
ethyl  bromide;  however,  pyrolysis  of  the  ethyl  bromide  caused  the  melt 
to  become  contaminated  with  carbon  particles.  Future  work  will  include  the 
growth  of  bromate-doped  KBr  from  melt  and  the  growth  of  pure  and  bromate- 
doped  KBr  from  solution. 

a.  3b  Optical  Absorption  of  KBr:  As  described  in  previous  reports,  the 
10.6  micron  bulk  absorption  coefficient  of  our  KC1  crystals,  grown  from 
purified  melts,  is  comparable  to  the  lowest  values  obtained  elsewhere1,2 
and  is  well  below  the  limit  of  reliability  for  our  present  calorimetric 
apparatus.  Therefore,  pending  improvements  in  that  apparatus  (see  section 
f  1)  we  are  currently  concentrating  on  measurements  in  pure  and  doped  KBr. 

As  described  in  the  last  report,  single  crystals  of  KBr  have  been  grown 
from  salt  which  was  purified  by  ion-exchange,  fractional  crystallization  and 
gas  treatment  of  the  molten  salt  with  HBr  and  Br2*  Specimens  for  calorimetry 
were  cleaved  from  these  crystals  without  subsequent  polishing  since  during 
our  previous  studies  of  KC1,  we  found  that  cleaved  samples  can  produce  better 
surfaces  than  our  present  polishing  technique. 

Initial  at.empts  to  measure  the  absorption  by  a  calorimetric  method  were 
not  successful  due  to  scattering  of  the  laser  radiation  to  the  thermocouple. 

This  scattering  was  caused  by  deterioration  of  the  crystal  surfaces  during 
the  time  the  samples  were  mounted  in  the  calorimeter.  The  results  of  sub¬ 
sequent  measurements,  presented  in  Table  1  and  Figure  1,  indicate  a  bulk 
absorption  (a)  of  0.0011cm"1  and  an  absorption  per  surface  (a)  of  0.00051cm"1. 

It  should  be  noted  that  positioning  the  sample  in  the  beam  so  as  to  eliminate 
scattering  was  very  difficult  and  that  the  uncertainty  of  this  set  of 
measurements  is  higher  than  for  previous  samples. 

For  comparison,  calorimetric  measurements  were  also  made  on  KC1  and  KBr 
cleaved  samples  purchased  from  Optivac.  These  are  used  as  seeds  for  crystal 
growth,  as  specimens  for  surface  preparation  studies,  and  as  standard  specimens 
for  comparison  between  our  results  and  those  of  other  laboratories. 

The  calorimetric  results  of  Optivac  crystals  in  the  "as-received"  condition 
(cleaved  at  the  factory)  are  listed  in  Table  2.  Figure  2  and  3  show  the  fit 
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of  the  data  to  the  equation 


82  =  01SL  + 


where  8  is  the  absorption  coefficient  per  centimeter,  i  is  path  length, 
a  is  the  absorption  per  surface  and  n  is  the  number  of  surfaces. 

The  results  for  as-received  KC1,  as  shown  in  Figure  2,  indicate  an 
“  °f  °-°0066cm-l  and  a  0  of  0.0021a,-!.  other  samples,  »ater  lapped  at  the 
factory,  had  such  cloudy  surfaces  that  they  were  not  measurable.  The 
same  was  true  of  the  water  lapped  KBr  specimens. 

For  comparison,  Table  3  lists  our  results  for  our  own  KC1  crystals, 
as  presented  in  the  previous  quarterly  report. 

The  measurements  on  Optivac  KBr  cleaved  samples  in  the  as-received 
condition  were  too  unreliable  -o  permit  extrapolation  to  estimate  a  and  a. 

This  was  because  of  excessive  scattering  of  radiation  (apparently  from  the 
surface)  which  made  reliable  measurements  difficult  to  make.  The  minimum 
surface  absorption,  from  Figure  3,  is  in  the  region  of  0.002cm^.  Evaluation 
of  KC1  and  KBr  samples  with  freshly  cleaved  surfaces  are  now  in  progress. 
Polishing^ techniques  are  being  re-evaluated  and  a  modification  of  Dyer's 
apparatus  may  soon  be  constructed.  Measurements  on  both  cleaved  and  polished 
KBr  crystals  will  be  made  during  the  next  quarter. 

— — °wth  KBr  Crystals;  It  has  been  suggested2  that  much  of  the 

10.6  micron  absorption  in  impure  KC1  and  KBr  crystals  might  be  due  to  the 
presence  of  chlorate  ions  in  KC1  and  bromate  ions  in  KBr.  Although,  it  is 
difficult  to  understand  how  bromate  ions  could  survive  in  molten  KBr  in  a 
HBr  atmosphere,  it  must  also  be  conceded  that  bromate  is  a  likely  impurity 
in  KBr  and  would  not  tend  to  be  removed  by  ion  exchange  or  fractional 
crystallization.  Accordingly,  since  bromate  is  a  strong  oxidizing  agent, 
it  wa  i  decided  that  a  more  reducing  reactive  gas  should  be  used  in  the  final 
stages  of  purification.  Of  the  gases  described  in  Lebl  and  Trnka's  original 
article  on  reactive  atmosphere  processing5  ethyl  bromide  seemed  most 
suitable.  Preliminary  experiments  with  this  reagent,  however,  have  been 
less  than  satisfactory  because  of  its  tendency  to  pyrolysis.  During  gas 
treatment  of  the  molten  salt,  a  heavy  black  deposit  of  carbon  “ormed  on  the 
silica  tube  of  the  crystal  puller  and  obscured  visibility  enough  to  seriously 
interfere  with  the  seeding  and  growth  process.  Moreover,  carbon  deposits 
were  also  apparent  on  the  inside  of  the  bubbling  tube,  indicating  that  the 
decomposition  of  C^Br  began  before  it  bubbled  through  the  melt.  A  few 
particles  of  carbon  were  found  in  the  resolidified  polycrystalline  salt; 
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Ue;e  Pr°bably  dlSl0d«d  tr»"  “»  tube.  There  are  „e  indications 

at  carbon  la  forced  In  the  molten  salt.  However,  since  cleaning  the  system 

between  the  gas  treatment  and  crystal  growth  stages  wonld  allow  possible 

contamination  of  the  salt  with  atmospheric  gases  and  residual  water  adhering 

to  the  meed  silica  liner,  It  seems  Inexpedient  to  continue  using  ethyl 

bromide.  Future  experiments  will  probably  Involve  reagents  such  as  CH  Br. 
and  CBr^  sublimed  into  a  carrier  gas.  2  2 

Work  is  also  progressing  on  the  growth  of  KBr  crystals  from  melts 
doped  with  KBr03  and  on  the  growth  of  KBr  crystals  from  aqueous  solution. 
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TABLE  1 


Crystal 
KBr  1  -  1 


Results  of  Absorption 

measurements 

of  KBr  at  10.6 

microns 

All  samples 

were  freshly 

cleaved 

Averagi 

Specimen  L(cm) 

8 (cm-1 ) 

BL 

A  .7359 

.00260 

.00204* 

.00199 

.00269 

.00212 

.00232 

.00182 

B  1.318 

.00203 

.00268 

.00242 

.00210 

.00277* 

.00192 

.00253 

.00181 

.00239 

.00164 

.00216 

C  . 4486 

.00  307 

.00138 

.00130 

.00296 

.00134 

.00294 

.00132* 

.00265 

.00119 

.00280 

.00126 

.00288 

.00129 

D  .2647 

.00548 

.00145 

.00141 

.00513 

.00136 

*Indicate  sample  repositioned  between  measurements. 
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TABLE  2 


Crystal 

OPT. 

KC1 


OPT. 

KBr 


Results  of 

Absorption 

measurements 

of  Optivac  KOI 

and  KRr* 

Specimen 

All  samples 

were  cleaved 

at  factory 

0L 

t* uu  I v  U  JL 

Al7Pl*Q  f> q 

L(cm) 

0 (cm-1 ) 

1  ('  a' 

«VCL  Clfc;“ 

axis) 

5.393 

.00194 

.01048 

.01135 

.00226 

.01217 

.00211 

.01140 

1  ('  b' 

axis) 

1.090 

.00464 

.00505 

.00509 

.00430 

.00/68 

.00507 

. 00553 

2  ('a' 

axis) 

2.040 

.00432 

.00881 

.00835 

.00453 

.00931 

.00339 

.00692 

2  ('b' 

axis) 

0.7117 

.00802 

.00571 

.00549 

.00737 

.00525 

.00655 

.00466 

.00890 

.00634 

1  ('a' 

axis) 

5.235 

.001207 

.006319 

.00556 

.000989 

.00518 

.000991 

.00519 

1  (V 

axis) 

1.114 

.00460 

.00513 

.00541 

.00549 

.00611 

.00447 

.00498 

2  ('a' 

axis) 

2.033 

.00258 

.00525 

.00517 

.00258 

.00525 

.00247 

.00502 

2  ('b' 

axis) 

.7407 

.00219 

.00216 

.00268 

.00393 

.00291 

.00329 

.00244 

r 


TABLE  3 


Summary  of  Results  of  Absorption  Measurements  at 
10.6  Microns  on  KC1  Crystal  3-7 


Sample 


Surface  Preparation 


u  (cm-1) _ o 


KC1  3-7  A,  B  and  C 
KC1  3-7  D  and  E 


Polished 
Freshly  cleaved 


0.00025 

<0.0001 


.0007 

'vO.OOOl 
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Figure  1  Plot  of  Absorption  Data  for  Freshly  Cleaved 
from  Crystal  KBr  1-1 


.6  .8  1.0  1.2  I 


LENGTH  ,  (cm) 
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LENGTH,  L (cm) 


LENGTH 


a. 3  Chemical  Vapor  Deposition  of  Bulk  GaAs  Crystals 
W.  R.  Wilcox,  E.  A.  Miller 

a. 3a  Introduction:  We  wish  to  summarize  the  effort  of  the  past  18  months 
which  was  directed  toward  growth  of  bulk  GaAs  by  the  open  tube  CVD  method, 
since  the  project  has  now  ended.  Particular  emphasis  will  be  on  the  last 
quarter's  work,  wherein  we  obtained  smooth  deposits  1  to  2  mm  thick  on  GaAs 

substrates . 

a. 3b  Experimental  Procedure:  An  open  tube  epitaxial  reactor  was  con¬ 
structed,  which  used  the  AsCl3-Ga-H2  deposition  system.  This  system  is  most 
often  used  in  the  semiconductor  process  industry  for  growth  of  high  purity  GaAs 
films  by  CVD  methods.  We  established  operating  conditions  favorable  for  the 
growth  of  thin  films,  but  found  that  these  conditions  were  not  suitable  for 
long  term  thick  layer  growth.  Next,  out  effort  consisted  of 

experimentally  determining  the  optimum  conditions  for  substrate  temperature, 

A?C13  mole  fraction  in  the  gas  substrate  orientation,  and  growth  time. 

a. 3c  Results  and  Discussion:  Several  long  runs  were  completed  which 
resulted  in  severe  substrate  etching.  This  was  accompanied  by  dendritic 
overgrowths  which  apparently  consumed  reactants  from  the  gas  stream  directly 
above  the  substrates.  By  reducing  the  AsC13  supersaturation  (lower  AsCl3  mole 
fraction  in  the  gas)  we  were  able  to  minimize  dendritic  growth.  It  seems,  however, 
that  given  enough  time  dendritic  overgrowth  occurs  regardless  and  acts  as  a 
factor  limiting  deposit  thickness.  For  example,  a  9  day  run  (CVD  #21)  gave 
smooth  deposits,  while  a  27  day  run  under  the  same  conditions  (CVD  #22)  resulted 
in  severe  etching  and  dendritic  overgrowths  (see  Table  I). 

From  results  of  previous  runs,  we  selected  operating  conditions  that  yielded 
deposits  1  to  2  mm  thick  in  7  days.  The  substrate  orientation  was  {100},  which 
reportedly  has  a  slower  growth  rate  than  other  orientations  .ested,  but  more 
uniform  deposits  usually  occur.  Figure  1  shows  the  results  of  this  run  (CVD 
#26);  note  the  clusters  of  GaAs  deposits  on  the  substrate  holder,  and  the 
blade-like  dendritic  growths  on  substrate  #1  located  at  the  tip  of  the  holder. 

Also. see  Figure  2  for  a  close-up  view  of  substrate  #1.  The  blades  on  substrate 
#2  (2nd  from  the  tip  of  holder)  are  shown  in  Figure  3. 

Substrate  #3  was  further  examined.  Note  the  appearance  of  large  individual 
grains,  Figure  1,  the  boundaries  of  which  were  revealed  by  inclusions  and  voids 
after  lapping  and  chemical  polishing.  The  orientation  of  these  grains  remained 
in  the  {100}  direction,  however.  All  deposits  in  this  run  had  many  small  voids 
and  inclusions,  possibly  due  to  small  air  leaks  in  the  system.  See  Figure  2 
for  example. 
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durlnsTa“e  1  S“™arlZe8  S°me  °f  th=  condition,  and  results  obtained 

during  this  program. 

Som  IR  absorptlon  anJ  „a 

data  were  taken  on  the  deposited  films.  Prom  run  CVD  #21.  „e  were  able  to 
obtain  a  sample  4  mils  thick  after  lapping  away  the  substrate.  Results  fro. 

1  h  a”P  ISble  Sh°“  ““  the  approXl.ately 

Of  he  same  guallty  as  typical  .it  grown  GaAs.  These  rather  poor  results 

ould  have  been  caused  by  the  rn.il  voids  and  Inclusions  observed  In  the  sample 
obT  rV  Sy8tOT  had  b“n  “<*  —  would  have  bT 

polycrystallinri^rT113^0113  ^  ^  °"  Sampl6S  Wlth  many  sma11  C-  0.5mm2) 
pclycrystalline  inclusions  in  them;  results  in  these  cases  wptp 

absorption  measurement  at  10  6u  a  eV6n  W°rSe'  An 

CVD  #26  Tb  b  T  °ne  °f  thP  bladeS  1"  nan 

^VU  ffZb.  The  absorbance  was  sn  hloh  (-hat  ,  ■> 

at  5  watts  ,  ,  ,  S  that  a  hole  ™  burned  through  the  sample 

at  5  watts  incident  power  Tn  ..aa- in  “ 

,  M  P  .  In  addition,  operation  at  the  lower  As Cl,  mole 

fractions  (necessary  to  produce  thick  films)  increases  the  tend  r 
incorporation  (mainly  Si)1  th  e  v.  tendency  for  impurity 

References  ’  ““  P-Pertles. 

1.  J.  v.  DiLorenzo,  J.  Crystal  Growth  17,  189  (19 72) 
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Figure  2  CVD  #26  Substrate  #1  4X  Figure  3  CVD  #26,  Stustrate  #2  iOX 


SUMMARY  OF  RUNS  3/1/73  -  11/30/73 
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b.l  Fabrication  of  Polycrystalline  IR  Window  Material 

S.  M.  Copley,  J.  M.  Whelan  and  V.  Rana 

The  aim  of  this  part  of  the  program  is  to  produce  samples  of  laser  window 
materials  in  polycrystalline  form  by  a  novel  hot  pressing  technique  involving  the 
use  of  volatile  sintering  aids.  This  technique  has  been  described  in  detail 
in  previous  quarterly  reports. 

Experiments  have  been  carried  out  by  hot  pressing  CdTe  powder  with  liquid 
Cd  as  a  sintering  aid.  CdTe  powder,  99.999  pet.  pure,  was  obtained  from  Research 
Organic/Inorganics.  This  powder  is  produced  by  coprecipitating  Cd  and  Te  from 
solution.  This  coprecipitation  is  carried  out  by  reducing  acidic  solutions  of 
Cd  and  Te  with  hydrazine  monohydrochloride  (t^H^HCl) .  The  resulting  precipitate 
is  washed  many  times  with  hot  water.  The  water  is  subsequently  driven  off  by 
heating  to  130°C  for  2  or  3  days  in  vacuum.  Average  grain  size  of  the  powder  used 
in  our  experiments  is  lym.  Extra  Cd,  to  be  used  as  a  sintering  aid,  added  to  the 
charge,  is  either  in  powder  form  (99.999  pet  pure,  -100  mesh)  or  as  shots  (99.999 
jet  pure) ;  both  obtained  from  Cominco  American  Incorporated. 

It  has  been  possible  to  prevent  any  loss  of  Cd  from  the  charge  by  using  the 
die  configuration  shown  in  Fig.  1.  By  utilizing  the  chill,  Cd  in  the  top  recess 
of  the  die  is  kept  at  a  slightly  higher  temperature  (AT  -  20°C)  than  the  Cd  in 
the  charge  so  that  it  provides  the  overpressure  of  Cd  required.  The  Cd  at  the 
top  is  prevented  from  escaping  by  covering  it  with  LiCl-KCl  eutectic  (melting 
point  355°). 

To  form  the  compact,  the  die  is  placed  in  the  furnace  inside  the  hot 
pressing  chamber  and  the  chamber  is  evacuated  to  a  pressure  of  30  micron  of 
Hg  or  better.  The  powder  is  then  hot  pressed;  two  procedures  of  hot  pressing 
have  been  employed.  In  the  first  procedure,  pressure  is  applied  to  the  die  plunger 
at  room  temperature  and  is  maintained  while  the  temperature  is  raised  at  a  rate  of 
approximately  2°C/min.  to  the  final  desired  temperature.  Curve  (a),  Fig.  2,  shows 
a  typical  curve  for  densification  obtained  by  this  procedure.  Microscopic 
examination  of  such  compacts  showed  that  Cd  was  not  able  to  disperse  throughout 
the  compact  uniformly,  but  was  segregated  at  several  places.  In  the  second 
procedure,  the  charge  was  first  heated  up  to  the  desired  temperature,  and  then 
pressure  was  applied.  All  the  compacts  obtained  in  this  manner  were  superior 
to  those  obtained  following  the  first  procedure,  curve  (b)  Fig.  2.  The  Cd  was 
much  better  dispersed,  and  after  wetting  the  grains  during  compaction,  the  extra 
Cd  was  squeezed  out  between  the  plunger  and  the  die  wall.  There  were  still,  however, 
certain  areas  of  compact  where  densification  was  poor.  Fig,,  3  shows  a  well  compacted 
area,  while  Fig.  4  shows  a  poorly  compacted  area  in  the  same  compact. 
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Severe  sticking  of  CdTe  compact  to  the  die  made  from  Speer  graphite  grade 
580,  was  observed.  There  was  also  evidence  of  Cd  effusing  through  the  walls  of 
the  die.  Using  Poco  graphite  grade  HPD-1,  these  problems  of  Cd  effusion  and 
sticking  have  been  almost  eliminated. 

CdTe  powder  used  in  our  experiments  has  been  found  to  contain  some  volatile 
matter.  Severe  outgassing  from  the  charge  was  noted  in  all  the  runs  (Fig.  5) . 

This  volatile  matter  could  be  any  of  the  reagents  used  in  synthesis  and  not 
completely  washed  away  (e.g.  hydrazine  monohydrochloride)  or  any  of  the  inter¬ 
mediate  compounds  not  completely  converted  to  the  final  product  CdTe  (e.g. 

Cd  (NH2)2  or  CdTe  N^)  .  Such  volatile  matter  if  not  removed  prior  to  compaction 
gives  rise  to  gases  that  can  get  trapped  inside  the  compact.  It  may  be  responsible 
for  the  voids  seen  in  the  poorly  compacted  samples.  It  may  also  prevent  the 
liquid  Cd  from  completely  wetting  all  the  particles.  It  is  also,  of  course, 
necessary  to  remove  this  volatile  matter  in  order  to  obtain  a  high  purity  compact. 

Efforts  are  underway  to  characterize  and  eliminate  the  volatile  impurities. 
Initially,  the  powder  will  be  heated  to  a  temperature  higher  than  the  boiling 
point  of  all  suspected  Impurities  and  its  effects  on  densif ication  seen.  Powder 
refinements  under  H2  gas  flow  will  also  be  tried. 


Relative  Density 


If  •  A  : 1*  >■  v'Hf-*"«y  ■>,.«  -! »$r*nr  ,1.  ,»  *,.■ 
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Fig.  2  (a)  1000  psi  applied  at  room  temperature  and  maintained,  while 

temperature  is  raised  at  the  rate  of  2°C/min  to  400°C 
charge  is  heated  to  400°C  and  kept  at  that  temperature  while 
3000  psi  is  applied  for  15  minutes. 


(b) 


*Based  on  the  theoretical  density  of  CdTe  of  5.86  gms/cm^ 
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Fig.  3  Fracture  surface  of 
to  400 °C 


a  compact  pressed  at  3000  psi  after  heating 


a  compact  pressed  at  3000  psi  after  heating 


C,;L  ^g£hanical  Behavior  of  III-V  and  II-VI  Compound 
S.  M.  Copley  and  V.  Swaminathan 


The  objective  of  this  investigation  is  to  characterize  the  mechanical 
behavior  of  various  candidate  window  materials  under  different  conditions 
and  to  identify  the  important  mechanisms  of  plastic  deformation  and  fracture. 

In  this  report  we  will  su-arize  briefly  the  results  to  date  and  indicate  the 
experiments  that  are  in  progress. 

c.la  Results  to  Date: 

—  ■■■■  H1Cr°  Structural  Characterizations.  We  have  employed  T.E.M.  to 
characterize  the  microstructure  of  several  types  of  GaAs  which  were  used  in 
the  investigation  of  mechanical  properties  and  also  several  types  of  high 
resistivity  GaAs  with  low  absoprtion  coefficients  investigated  by  Steier  and 

co-workers  as  part  of  this  program.  The  results  of  this  investigation  were 
reported  in  Quarterly  Technical  Report  H. 

CrlaJ_Momalo.us  Diffraction  Effects;*  Diffuse  intensity  streaks  and  extra 
spots  have  been  observed  in  transmission  electron  diffraction  patterns  of  gallium 
arsenide  single  crystals.  The  configuration  of  these  diffuse  streaks  and  extra 
spots  and  their  intensity  have  been  studied  as  a  function  of  crystal  orientation 
an  temperature.  The  geometry  of  the  diffuse  streaks  and  extra  spots  can  be 
interpreted,  respectively,  in  terms  of  the  Intersections  of  (110)  sheets  and 

SPlk<iS  °f  dlffuSe  lntenslty  1"  reciprocal  space  with  Ewald's  sphere.  The 
presence  of  {110}  sheets  of  diffuse  intensity  is  explained  on  the  basis  of  thermal 
use  scattering  of  electrons  by  low  frequency  acoustic  phonons  with  wave 
vectors  lying  on  (110)  planes.  The  presence  of  <111>  intensity  spikes  is  also 
attributed  to  thermal  diffuse  scattering  and  is  explained  as  a  combined  effect 
Of  lattice  symmetry  and  phonon  anisotropy.  This  work  was  presented  at  the 
Third  Conference  on  High  Power  Infrared  baser  Window  Materials,  and  will  be 
published  in  the  conference  proceedings. 

f  ~~  Stress-Straln  Savior:  Stress-strain  curves  have  been  obtained 

or  i-doped  GaAs  single  crystals  to  <100>  and  <111>  at  temperatures  ranging 
rom  250  to  550’G.  I„  both  orientations,  slip  is  shown  to  occur  only  on  {111} 

<  >.  Repeated  yielding  experiments  on  GaAs  crystals,  both  undoped  and  doped 

”  '  ’  ^  ^  Z"  lndlCSte  that  S1-<‘°Ped  <^tals  (n-type)  have  a  greater  yield 


*W„rk  described  in  this  section  was  carried  out  in  collaboration  with  Professor 
G.  H.  Narayanan,  Dept,  of  Materials  Science,  U.S.C. 
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stress  than  undoped  crystals,  and  undcped  crystals  have  a  greater  yield  stress 
than  Zn-doped  crystals  (p-type) ;  Cr-dcping  does  not  significantly  affect  the 
yield  stress.  The  effect  of  impurities  is  attributed  to  dynamic  hardening 
through  the  formation  of  tetragonal  defects  involving  vacancy -solute  pairs. 
Extensive  crack  formation  due  to  deformation  is  not  observed.  Compression 
specimens  normally  failed  by  breaking  into  several  columnar  fragments  parallel 
to  the  stress  axis.  The  results  of  this  work  were  presented  at  the  Third 

High  Power  Infrared  Laser  Window  Materials  Conference  and  will  be  published  in 
the  Conference  proceedings. 
c.lb  Work  in  Progress: 

— -lb ' 1 — I™Purlty  Hardening:  Both  theoretical  and  experimental  work  is 
currently  in  progress  with  the  £oal  of  relating  the  observed  impurity 
hardening  to  the  presence  of  tetragonal  vacancy-solute  pairs.  GaAs  single 
crystals  doped  with  Si  in  the  range  of  1018  -  1019/cm3  will  be  given  different 


annealing  treatments  prior  to  deformation.  These  annealing  experiments  will  be 
based  on  results  by  Kung  and  Spitzer1.  Samples  will  be  annealed  at  1100°C  for 


15  minutes  and  they  will  be  aged  at  700°C,  600°C,  and  450°C  for  specified 
periods  of  time.  The  effect  of  different  concentrations  of  such  solute-vacancy 
pairs  on  mechanical  behavior  will  be  studied.  These  experiments  will  be  repeated 
with  Te  doped  GaAs.  Additional  experiments  on  Si-doped  with  P  and  B  are  also 
being  carried  out  to  determine  whether  the  solute-vacancy  pair  mechanism  can 
explain  the  impurity  effect  in  the  elemental  semiconductors. 

Annealing  experiments  will  also  be  carried  out  on  GaAs  containing  Si  >  1019/cm3, 
Recent  work  at  U.S.C.  by  Narayanan  and  Kachare  has  shown  that  annealing  treatments 


on  such  highly  doped  crystals  can  cause  the  formation  of  a  high  density  of 

2 


prismatic  dislocation  loops  of  the  vacancy  type2.  The  annealing  specimens  will 
be  tested  in  compression  so  that  the  effect  of  such  prismatic  dislocation  loops 
on  strength  can  be  determined. 

A  special  loading  ram  is  being  designed  to  carry  out  compression 
experiments  at  room  temperature  with  precision  alignment.  The  possibility  of 
twinning  as  a  low  temperature  deformation  mechanism  is  being  investigated. 

cj.b.2  Cold  Work  and  Process  Annealing:  We  have  established  conditions 
for  producing  a  recrystallized  substructure  in  GaAs  single  crystals  (see 
Quarterly  Technical  Report  H) .  It  is  believed  that  recrystallized  material 
may  well  show  improvements  in  absorption  in  comparison  to  melt  grown  single 
crystals.  This  is  because  the  recrystallization  process  may  lead  to  a 
redistribution  of  impurities  along  grain  boundaries.  Recrystallization 
experiments  will  be  carried  out  on  high  resistivity  Cr-doped  GaAs  obtained 
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irom  Bell  and  Howell.  Specimens  will  be  deformed  to  20%  at  550°C  and  then  they 
will  be  annealed  at  110"C  for  10  minutes.  This  treatment  will  produce  the 
recrystallized  structure.  Absorption  as  a  function  of  wavelength  will  be 
measured  before  and  after  recrystallization. 

cab. 3  Mechanical  Behavior  of  CdTe :  During  the  renewal  period  of  this 
project  a  basic  investigation  of  the  mechanical  behavior  of  CdTe  will  be  carried 
out.  The  material  will  be  obtained  from  II-VI,  I„c.  Our  main  approach  will 
be  to  emphasize  the  effect  of  thermal  and  thermo  mechanical  treatments  on  the 
defect  structure  and  absorption  of  CdTe.  The  structure  of  specimens  annealed 
under  specified  partial  pressures  of  Cd  and  Te  will  be  studied  in  collaboration 
with  Dr.  Kroger.  The  Investigation  of  the  dlfferenct  microstructures  will  be 
carried  out  by  transmission  electron  microscopy.  Thin  foils  of  CdTe  have  already 
been  successfully  prepared.  The  orientation  and  temperature  dependence  of 

plastic  deformation  will  be  studied.  Mechanisms  of  deformation,  especially 
twinning,  will  be  investigated. 
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d.l  Dielectric  Constant  Measurements 


Clarence  R.  Crowell,  S.  Joshi 

Effort  during  this  quarter  was  directed  towards  preparing  more  GaAs 
samples  and  confirming  the  dielectric  constant  measurements  performed  on  the  first 
sample  and  reported  in  September,  1973. 

A  second  sample  was  prepared  on  another  wafer  of  the  same  material.  More 
precise  marking  for  defining  the  active  area  was  adopted.  The  sample  was 
prepared  by  the  same  method  as  reported  previously.  Capacitance  vs.  voltage 
characteristics  and  capacitance  vs.  frequency  characteristics  were  measured  and 
were  found  to  be  the  same  as  for  the  first  sample.  The  variation  of  capacitance 
vs.  frequency  at  liquid  nitrogen  temperature  reported  earlier  was  eliminated  by 
better  probe  contacting  technique  and  a  flat  capacitance  vs.  frequency  curve 
was  obtained  at  that  temperature. 

The  measurements  of  thickness  and  area  on  the  sample  gave  an  effective 
dielectric  constant  about  10%  higher  than  reported  earlier.  The  thickness  of 
the  first  sample  was  remeasured  using  present  techniques  and  was  found  to  be 
in  agreement  with  a  revised  dielectric  constant  of  approximately  13.6. 

A  need  was  felt  to  confirm  this  result  on  more  samples  and  in  particular 
to  check  the  uniformity  of  a  given  sample.  Accordingly,  a  sample  was  prepared 
from  WA-100  GaAs  material  and  gold  electrodes  were  deposited  to  define  dielectric 
capacitors.  Area  and  thickness  measurements  on  these  capacitors  gave  an  effective 
dielectric  constant  which  varies  about  3%  from  capacitor  to  capacitor. 

From  these  results  it  is  felt  that  the  sample  flctness  and  thickness 
measuring  techniques  still  need  to  be  improved  to  obtain  the  dielectric 
constant  with  better  accuracy. 

Preparation  of  dielectric  capacitors  on  CdTe  material  for  measuring  the 
dielectric  constant  is  also  underway. 


High  temperature  Hall  effect  measurements  on  lightly  doped  CdTe  are 
being  repeated  to  check  on  a  small  discrepancy  between  our  data  and  those  of 
Smith.  The  results  will  be  compared  with  those  of  similar  measurements  in 
cooled  crystals  in  order  to  see  what  happens  to  the  native  donors  upon  cooling. 

We  intend  to  grow  variously  doped  crystals  for  local  mode  studies  by 
Dr.  Spitzer  and  resistivity  Hall  effect  and  10 p  absorption  for  ourselves. 

A  model  to  explain  the  properties  of  Ag-doped  CdS  has  been  developed, 
and  the  observed  Hall  and  resistivity  data  are  being  interpreted  on  the  basis 
of  this  model. 
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L.  DeShazer  and  M.  Bass 


During  this  quarter  we  began  to  assemble  the  pulsed  CO  TEA  laser 
for  use  in  the  10.6pm  high  peak  damage  studies.  A  few  of  the  pieces  have 
been  fabricated  and  a  few  components  ordered.  Most  of  the  fabrication, 
assembly  and  testing  will  be  completed  during  the  next  quarter. 


e . 1  Theoretical  Studies  of  Absorption  Mechanisms  in  IR  Window  Materials 


R.  W.  Hellwarth  and  M.  Mangir 

e.l.a  Determination  of  the  Mechanism  of  Multi-Phonon  Absorption  in 
Pure  Crystals 

This  quarter  was  devoted  mainly  to  the  search  for  a  method  by  which 
one  can  determine  unambiguously  from  experiment  the  relative  contribution  of 
multi  phc.ion  absorption  in  a  given  material  of  the  two  possible  mechanisms: 
a)  anharmonic  interionic  forces  and  b)  the  nonlinear  dependence  of  the  electric 
dipole  moment  on  ionic  displacements,  i.e.,  the  "nonlinear  moment".  We  believe 
we  have  discovered  the  first  practical  method  to  achieve  this.  It  is  described 
in  detail  in  Appendix  I  which  is  to  appear  in  the  Proceedings  of  the  Third 
Conference  on  High  Power  Laser  Window  Materials  (AV  *Y>rce  Cambridge  Research 
Labs,  1973)  ea.  by  Carl  Pitha.  As  is  pointed  out  there,  so  far  the  necessary 
data  for  our  method  exists  only  for  LiF  (from  which  we  concluded  that  mechanism 
(a)  dominates  in  pure  LiF  at  room  temperature) .  We  are  currently  trying  to 
collaborate  with  experimenters  at  USC  and  elsewhere  in  making  the  measurements 
necessary  to  determine  the  multi-phonon  absorption  mechanisms  in  III-V  and 
II-VI  window  materials.  For  these  materials  the  mechanism  is  especially  unclear 
at  this  point. 


Appendix  for  e.l 


A  METHOD  FOR  MEASURING  THE 
CONTRIBUTION  OF  THE  NONLINEAR  DIPOLE. 
MOMENT  TO  MULTI -PHONON  ABSORPTION" 


ROBERT  HELLWARTH  and  MET IN  MANGIR 
Electrical  Engineering  and  Physics 
University  of  Southern  Calif. 
Los  Angeles,  Calif.  90007 


Abstract 


A  method  is  proposed  by  which  one  may  determine  the  separate 
contributions  of  the  anharmonic  interionic  forces  and  of  the  nonlinear 
dipole  moment  of  any  material  to  the  integral  over  IR  frequencies  of 
the  imaginary  part  of  its  dielectric  function  (times  the  frequency). 
With  this  we  show  from  existing  data  that  the  contribution  of  multi¬ 
phonon  absorption  to  this  integral  for  LiF  is  predominantly  from 
anharmonic  forces. 


1.  INTRODUCTION 


In  order  to  calculate  the  multi-phonon  absorption  in  a  given 
material,  one  must  first  know  the  relative  contribution  of  the  two 
possible  mechanisms:  A)  anharmonic  interionic  forces,  and  B)  the 
nonlinear  dependence  of  the  electric  dipole  moment  on  ionic  displace¬ 
ments  (or  the  "nonlinear  moment").  For  alkali-halides ,  it  has  been 
found  consistent  in,  albeit  approximate,  calculations  to  neglect 
mechanism  (B) .  However,  for  some  promising  high  power  laser  window 
materials,  such  as  ZnSe,  Cdfe,  and  GaAs,  the  relative  importance  of 
the  two  mechanisms  remains  uncertain. ^ 
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Wc  propose  here;  an  exper  i  menta  1  method  for  assessing  these  mechanisms 
in  any  material.  The  method  is  based  on  the  well  known  fact  that 
the  frequency  shift (s)  and  width (s)  of  the  fundamental  infrared 

lattice  absorption  resonance  (s)  are  insensitive  to  tire  nonlinear 

/" 

moment,  whereas,  multi-phonon  difference-  and  summation-absorption 
bands  are  sensitive  to  both  mechanisms. ^  At  the  present  the  only 
material  for  which  enough  data  exists  to  apply  the  method  is  LiF, 


and  we  find  that  the  anharmonic  force  mechanism  (A)  is  indeed  mainly 
responsible  for  the  more  prominant  multi-phonon  features  in  the 
absorption  spectrum  cf  LiF.  (See  Section  4) . 

Our  method  requires  the  measurement  of  the  imaginary  part 
01  the  dielectric  constant,  throughout  the  infrared,  and  at  two  or 
more  temperatures.  From  this  measurement  one  obtains  an  experimental 
determination  of  the  integral  in  the  following  exact  sum  rule,  in  esu 
(proved  in  Appendix  A) 

2 

/"  dWWE2(U)  -  raj.  <(3Mx/3rici)2>/ina  sec’2  ,  CD 


in  which  M  is  the  operator  representing  the  electric  dipole  moment 

X 

along  the  principle  dielectric  x-axis  for  which  was  measured. 

For  optically  isotropic  materials  this  x-axis  can  be  any  axis  and 

all  principle  dielectric  constants  are  equal.  The  r^a  are  the 

spatial  coordinates  (i=x,y,z)  of  all  the  N  ions  (a=l,  2,...N),  with 

masses  m  ,  inside  a  macroscopic  volume  V  of  the  crystal.  The 
a 

brackets  <>  indicate  the  quantum  average. 


The  usefulness  of  (1)  lies  in  the  fact  that  the  contrioution 


of  the 
can  be 
To  see 


averaged  expression  of  the  linear  term  in  the 

easily  calculated  exactly  for  many  important 

this,  we  first  write  each  ion  coordinate  as 

r  eR  +x. 
ict  ia  ’ 


dipole  moment 
types  of  crystals. 

(2) 
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the  sum  of  its  average  value  R-  E  <r.  >  (for  some  temperature)  and 

1(1  la 

its  deviation  x.  from  this  average  value.  Then  the  dipole  moment 

in  a  volume  V  of  the  crystal  at  this  temperature  can  be  expanded  in 

a  Taylor  series  in  the  x.  : 

ia 

Mx  "  Eiaexiaxiat,s!:iajBfxiajBxiaxjB+---  ’  t3) 

where  the  e  •  and  f  .  .  „  are  the  Taylor  series  expansion  coefficients, 

xia  xiaj  6 

the  former  being  an  "effective  charge"  tensor  for  the  at^1  ion.  The 
terms  in  M  third  and  higher  order  in  tr.e  x-  may  be  assumed  to 
contribute  less  to  (1)  than  the  quadratic  term  whenever  the  quadratic 
term  gives  rise  to  sum-and  difference-absorption  bands  that  contri¬ 
bute  less  to  the  integral  in  (1)  than  does  the  fundamental  absorption 
line(s)  arising  from  the  leading,  linear  dipole  moment,  term.  (In 
fact  for  all  potential  IR  window  materials  we  may  neglect  higher  than 
quadratic  terms  in  M  when  computing  the  RHS  of  (1)  to  within  the 
experimental  errors  in  determining  the  LHS)  . 


2.  CRYSTALS  WITH  LINEAR  DIPOLE  MOMENTS 


If  the  mechanism  of  anharmonic  forces  dominates  that  of  the 

nonlinear  moment,  we  may  neglect  all  but  the  term  linear  in  x^a  in 

the  expression  (3)  for  the  dipole  moment.  Then  we  need  study  only 

the  e  •  ,  whose  forms  depend  on  the  symmetry  of  the  site  at  R.  . 

xia  '  ia 

For  ions  at  sites  having  point  symmetries  T,  T,  ,  T , ,  0,  or  0  (such 

n  a  h 

as  the  ions  in  the  alkali  halides,  GaAs,  ZnSe  and  CdTe)  , 

xia  a  xi  ' 

where  e  is  called  the  "effective  charge"  of  the  ion  a  and  6  .  is 
a  °  xi 

the  Kronecker  delta.  For  diatomic  cubic  crystals  such  as  the  afore¬ 
mentioned,  the  two  effective  charges  must  be  related  by 

* 

el  =  'e2  =  e  (5) 
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in  0rdCr  U,at  "°  m0,"e,U  from  (3)  when  all  urns  a,c 

uniformly  translated. 

If  all  the  prominent  features  of  c2Cm)  arise  from  a  linear 
dipoie  moment ,  acting  with  enharmonic  interionic  forces,  in  aUali- 

haHde  ^  °tHer  -»ic  crystals,  then  we  see  from  fj,.C5j 

that  (1)  reduces  to 


CO 

/Q  ->■  2  7i^  pe  y  , 


(6) 


where  p  is  the  number  density  of  molecules  and  -l+m  -1  is  the 

a  B 


molecular  reduced  mass.  In  this  case,  furthermore,  t“he  crystal 

must  exhibit  a  dielectric  function  e2  that  obeys  (6)  at  different 
t  crop  c  r  ci  t  u  r  g  s  cind  dpns  i  t  i  p  q 

,  ensures,  with  an  essentially  invariant  effective 

ge  e  .  This  can  be  seen  as  follows.  As  the  temperature  is 
raised  and  the  crystal  expands,  the  Taylor  series  (3)  for  the  dipole 


moment  must  be  re-expressed  in  Dowers  nf  the  a 

powers  of  the  deviations  x.  '  about 

la 


the  new  average  ionic  positions  R.  »  =  R  +,  „  q  ,  . 

la  ia6Kia*  Substituting 


x  =  x,  ’  +6  R 


If* 


la 


ia 


(7) 


in  (3)  gives  for  the  "effective  charge"  coefficients  of  the  linear 
term  in  the  dipole  moment  expansion,  at  the  new  temperature, 


xi  a 


=  e  .  +  Z  .  f 

*v  -i  rv  -I  /-»  A 


xia  jB  xiajB6JB 


•  d6  R-  +. . . 


(8) 


But  6RjfJ  is  of  the  order  of,  or  smaller  than,  the  rms  fluctuations 

in  v  TU _ r  ..  .  ^ 


,  X.uo  iiutiuaiions 

in  xia.  Therefore,  if  the  nonlinear  moment  term  can  be  neglected  in 

f  1  tllon  m  j.1.  ■» 


- - ^  ^ 

(3),  then  so  can  the  change  in  effective  charge,  proportional  to 
(SR,  in  (8)  .  QED. 


3.  CRYSTALS  WITH  NONLINEAR  DIPOLE  MOMENTS 


If  the  imaginary  part  of  the  dielectric  constant  fails  to  obey 
(6),  then  it  can  only  be  because  the  nonlinear  quadratic,  or  "f" 
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term  in  tnc  dipole  moment  (3)  must  be  contributing  to  the  average  in 
(J).  To  see  what  a  deviation  from  (6)  means,  one  must  construct 
a  theory  for  the  form  of  the  f  -  coefficients,  and  use  experimental 


values  of  the  LHS  of  (1)  to  determine  the  parameters  in  the  theory. 

2 


Fortunately,  the  averages  over  the  (fx)  terms  required  in  (1)  can 
always  in  principle  be  performed  because,  to  the  desired  order  of 
accuracy,  the  lattice  may  be  assumed  to  be  perfectly  harmonic.  We 
illustrate  the  foregoing  procedure  by  using  the  independent -molecule 
model  of  the  lattice  that  has  been  used  extensively  in  estimating 
multi-phonon  absorption . C In  this  model,  each  one-dimensional 
molecule  in  the  crystal  contributes  independently  to  the  total 
dipole  an  amount 

,  ,  2 

mx  -  ex  +  ijfx  (9) 

where  x  is  the  deviation  from  the  equilibrium  interatomic  distance 
at  the  temperature  under  consideration.  Equation  (1)  then  reduces 
to 


!0  dwme2  (w)-*2iT2p(e2  +  f2<x2>)/y 


(10) 


At  another  temperature  (indicated  by  primes)  the  frequency  integral 
becomes 

2  TT2  p  ’  [(e  +  f6R)2  +  f2<x’2>]/p  .  (11) 

The  linear  expansion  <5R  and  densities  p  and  p'  are  known,  and  <x2> 
is  easily  calculated  for  the  harmonic  oscillator  at  the  two  tem¬ 
peratures,  assuming  its  mass  to  be  y  and  its  frequency  to  equal 
“TO*  There  fore  the  two  expressions  (10)  and  (11)  can  be  solved 
for  the  two  unknowns  e  and  f.  Then  the  contribution  to  the  integral 
of  multi-phonon  absorption  arising  from  the  nonlinear  moment  is 
simply  the  value  of  the  second  term  in  (10)  or  (11).  Subtracting 
this  value  from  the  total  multi-phonon  contribution  measured  gives 
the  anharmonici ty  contribution  to  the  integral  in  (10)  or  (11). 
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In  this  manner  the  multi -phonon  absorption  mechanisms  can  be 
separated  with  fair  accuracy.  Introducing  more  sophisticated 
forms  for  Mx  than  (9)  into  more  realistic  lattice  models,  and  using 
the  lattice  shell  models"  that  have  been  developed  to  fit  neutron 
scattering  data,  should  give  extremely  accurate  separation  of 
mechanisms  once  accurate  data  on  e2  (to)  becomes  available.  However, 
we  believe  that  even  the  simple  independent -molecule  model,  used 
judiciously  in  (1),  can  yield  a  very  good  picture  of  the  relative 
importance  of  the  two  mechanisms.  We  now  illustrate  its  use  in  the 
following  section. 

4.  APPLICATION  TO  LITHIUM  FLOURIDE 


The  only  crystal  for  which  ^(w)  has  been  measured  over  a 
wide  frequency  range  at  several  temperatures  is  LiF.  M.Klier4  has 
measured  its  by  transmission  measurements  in  the  ranges  8  to 
14pm  and  40  to  400pm,  at  300,  22,  and  -196°C. 

Kachare,  et  al.,5  have  analyzed  the  reflection  data  of 
6 

Jasperse,  et  al.,  by  classical  dispersion  analysis  to  give  values 
for  e2  in  the  range  omitted  by  Kleir  and  at  temperatures  7.5,  85, 
295,  420,  and  600°K.  Kachare  (private  communication)  later  veri¬ 
fied  these  results  by  a  Kramers- Kronig  analysis.  We  have  interpo¬ 
lated  this  data  to  combine  it  with  Klier's  and  evaluate  the 
integral  of  (1)  at  77,  295,  and  573°K. 

The  results,  with  number  densities  from  ref.  7,  are  given  in 


Table  1. 


Table  I 


Tempe 


rature  (°K) 

10"2  8/oduujC2  (uO  sec’" 

10_22p  cm 

77 

3.65 

6.20 

295 

3.64 

6.12 

573 

3.72 

5.92 
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We  esti,nate  the  experimental  accuracy  of  the  integral  to  be  ±K 
from  the  errors  estimated  in  refs.  4  and  5.  If  this  estimate  is 
correct,  the  experimental  results  are  not  quite  consistent  with  a 
pure  linear  dipole  moment,  once  the  density  changes  are  accounted 
for.  If  we  solve  (11), taking  u  to  be  8.44  x  icr24gm  (the  reduced 
mass  of  Li  and  F) ,  taking  the  molecular  vibration  frequency  to  be 
the  observed  fundamental  lattice  absorption  frequency  (given  in 
r_f.  5  as  317.5,  307.5,  and  300  cm  1  respectively  at  the  temperatures 
in  Table  I),  and  assuming  the  values  of  the  integrals  in  Table  I 
to  be  exact,  we  find  that  the  f2<x^>  term  contributes  2.21  of  the 
integral  at  573°K  whereas  10  ±  1%  of  the  integral  is  observed  to 
arise  from  multi-phonon  absorption.  If  the  measured  values  for 
the  integrals  in  Table  I  were  in  error  so  as  to  underestimate  the 
nonlinear  moment  effects,  we  estimate  that  the  f2<x2>  term  could 
contribute  as  much  as  a  third  of  the  multi-phonon  part  at  573°K. 

The  values  for  295°K  suggest  that  the  nonlinear  moment  is  relatively 
less  important  at  the  lower  temperatures  where  the  multiphonon 
contributions  are  less  altogether .  (From  (ID)  e  =  1.04  electron  charge.) 

In  summary,  the  uncertainties  in  the  data  and  in  the  approx¬ 
imations  of  (10)  and  (11)  probably  allow  us  to  conclude  no  more 
than  that  the  nonlinear  moment  gives  a  minor  if  not  negligible 
contribution  to  the  multi-phonon  absorption. in  LiF  as  it  is  mani¬ 
fested  in  the  integral  in  (1).  We  hope  that  similar  measurements 
will  be  made  in  materials  like  ZnSe,  CdTe  and  GaAs  for  which  the 
mechanisms  are  so  much  more  in  doubt.  The  application  of  the  exact 
relation  (1),  or  approximate  forms  like  (10)  and  (11),  to  such  data 
should  be  able  to  establish  the  relative  importance  of  anharmonicities 
and  nonlinear  moments  to  multi-phonon  absorption  in  these  materials. 
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APPENDIX  I 


According  to  a  well-known  quantum  mechanical  formula  (in  esu)^ 
e2M  =■  ^2VlEm(nOVl.n)|<m|Mx|n>|26(fl„.En+Em)  (A1) 

where  P^  is  the  probability  of  the  material  inside  V  being  in  state 
|m>  with  energy  Integrating  (Al)  by  wdw  from  to « replaces 

the  6-function  in  (Al)  by  (E^EJ/di.  Then  the  Pn  terms  give  the 
same  result  as  do  the  Pm,  and  may  be  dropped  if  the  latter  are 
doubled.  Next,  use  the  identity 

VEm'En)  I  < m  |  Mx  |  n >|2  =  i5 < m  |  [  [ H , ]  ,  Mx]|m>  (A2) 

and  then  evaluate  the  commutator  using 

H  =  hi  p  2/m  +U  fA3l 

3a  2a  a  • 

Since  both  the  full  lattice  potential  energy  U  and  M  are  functions 
of  the  ion  positions  only,  U  does  not  contribute  to  the  commutator. 
Using  the  commutator  identity 

. .  (A4) 


M  ,p  .  =  rR3iM  /3x. 

x  ia  x  i  a 


in  (A3)  leads  then  directly  to  the  desired  sum  rule,  equation  (1)  . 
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f-1  aSaUgg!i^^  Measurements  on  .... 

—  H‘  Si:eler>  c-  P.  Christensen.  R.  Joiner 

In  the  preceding  quarter  we  reported  irrevereible  increases  in  the 

absorption  of  high-resistivity  GaAs  produced  by  heating  the  material  i„  air  to 

a  temperature  of  Wc.  At  that  time  it  was  evident  that  both  surface  and  bulk 

effects  contributed  to  the  observed  change  in  room-temperature  absorption. 

Attempts  to  determine  the  thickness  of  the  highly-absorbing  surface  !ayer  by 

use  of  controlled  chemical  polishing  or  etching  have  been  unsuccessful  due  to 

formation  of  an  opaque  surface  film  during  the  polishing  or  etching  process  when 

less  than  5u  -  lOp  of  material  is  removed.  It  has  been  found,  however,  that  a 

cleaning  process  using  the  organic  solvents  trichlorethylene,  acetone,  and  methyl 

alcohol,  followed  by  wiping  the  surfaces  with  a  soft  sw*h  t  u 

wicn  a  soft  swab  or  cloth  is  sufficient 

to  remove  a  siaable  fraction  of  the  absorbing  surface  layer. 

The  absorption  mechanism  in  the  surface  laver  is  nnr  n 

ayer  is  not  well  understood.  It 

mental  vibrations  with  frequencies  as  high  as  1200  »-««.«  responsible  for 

Cbe  observed  loss,  however,  an  attempt  to  measure  surface  oxygen  concentration 

by  electro-chemical  techniques  with  the  assistance  of  Professor  J.  M.  Whelan 
was  inconclusive. 

Bulk  absorption  in  samples  subjected  to  several  thermal  cycles  has  been 

found  to  increase  after  each  heatine  Rpli  T  u 

ing.  Bell  Laboratories  sample  55-1A,  described 

in  the  previous  report,  showed  an  increase  in  bulk  absorption  from  .006  cm"!  to 

•018  cm-1  after  four  successive  400°C  heating  cycles.  However,  the  increased 

bulk  absorption  was  not  accompanied  by  significant  changes  in  electrical 

properties.  Resistivity  and  Hall  measurements  made  after  the  heating  cycles 
yielded  p  -  2.8  x  lo*  n  -c„  and  ^  ,  3350  cn,2/c_scc.  ComparIson  ^  ^ 

LaI>8  datS  °n  SlmUar  SampUs  #*)  shows  a  possible  factor  of 

-  reduction  in  mobility  to  be  the  only  effect  of  heating  as  far  as  electrical 

parameters  are  concerned.  The  irreversible  changes  caused  by  heating  are  thus 


still  a  mystery. 

The  temperature  dependence  of  multiphonon  absorption  in  GaAs  near  lOy  has 
been  calculated  using  the  formulas  of  McGill  et  al.  and  Bendow^.  Using 

the  experimental  results  of  Chang  et.al.^and  Horrigan  and  Deutsch^5'*  to  obtain 
the  necessary  constants,  we  find  good  agreement  between  experiment  (Q.R.  #5)  and 
theory  for  temperatures  below  280°C.  However,  it  is  clear  that  other  absorption 
mechanisms  could  also  exhibit  the  observed  weak  temperature  dependence  in  the 
temperature  range. 

Absorption  spectra  for  two  additional  samples  have  been  obtained  in  a 
continuation  of  our  GaAs  evaluation  program.  Sample  44-3A,  a  carbon-doped 
crystal  obtained  from  Bell  Laboratories  (see  Q.R.  #4  for  resistivity  and 
mobility),  exhibited  an  absorption  spectrum  nearly  identical  to  the  other  very 
low-loss  samples  that  have  been  measured  (55-1A,  WA-1000,  WORM-1).  A  second 
high  resistivity  sample,  doped  with  Si  and  Cr,  was  obtained  from  Hoboken  Lab¬ 
oratories  of  Belgium  and  showed  an  absorption  spectrum  very  similar  to  sample 
EMC  6050T  reported  pre\iously. 

We  note  that  the  four  lowest-loss  samples  in  the  evaluation  program 
representing  widely  different  growth  and  comparison  techniques,  exhibit  nearly 
identical  absorption  spectra.  This  suggests  that  absorption  in  the  best 
currently  available  GaAs  may  be  limited  by  an  intrinsic  loss  mechanism.  However, 
the  observed  irreversible  heating  effects  require  explanation  if  this  hypothesis 
is  to  hold.  Experiments  are  now  n  preparation  which  will  allow  annealing  of 
samples  in  a  controlled  atmosphere.  These  experiments  should  provide  further 
information  on  the  origin  of  the  irreversible  effects. 
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j-ii_A1kali  Halide  Surface  Studies  with  Acoustic  Probe  Techniques 
J.  H.  Parks,  D.  A.  Rockwell,  T.  Colbert 

During  this  quarter  we  successfully  demonstrated  that  a 
fluid  coupling  process  ^  could  be  used  to  transfer  surface 
acoustic  waves  from  crystalline  quartz  onto  non-piezoelectric 
KC1.  The  establishment  and  control  of  this  technology  allowed 
us  to  make  the  first  measurements  of  10. 6y  absorption  in  KC1 
samples  using  the  techniques  of  acoustic  temperature  detection. 

f.2a  Fluid  Coupling  of  Acoustic  Surface  Waves 

initial  fluid  coupling  experiments  indicated  that  several 
liquids  could  be  useful.  These  included  water,  xylene,  and  a 
eutectic  Ga : In  binary  alloy.  Xylene  was  chosen  as  the  coupling 
fluid  for  KC1  on  the  basis  of  minimum  phase  drift  and  coupling 
stability.  The  coupling  technique  is  illustrated  in  Figure  1. 

The  transmitting  and  receiving  transducers  were  constructed  on 
a  single  piece  of  quartz  to  insure  proper  orientation  and 
positioning  with  respect  to  each  other.  The  KC1  samples  were 
obtained  from  Janos  Optical  Corp.  and  had  a  single  surface 
optically  polished  to  a 1 low  propagation  of  the  acoustic  surface 
wave.  An  RF  signal  of  120  MHz  and  500  mV  rms  was  applied  to  the 
transmitting  trans^'cer  and  ~1  mV  rms  signal  was  detected  at 
the  receiving  t1  cer .  This  low  coupling  efficiency  is  due 

to  the  acoustic  impedance  mismatch  at  the  fluid  coupling  inter¬ 
faces.  This  could  be  increased  somewhat  if  we  replaced 
quartz  with  another  piezoelectric  material  having  acoustic  pro¬ 
perties  closer  to  those  of  KC1,  such  as  Bismuth-Germanum-Oxide 
(BGO) .  The  important  point,  however,  is  that  the  coupling 
technique  doe.  perform  adequately  and  the  received  signal  is 
of  sufficient  strength  to  allow  accurate  determination  of  the 
surface  wave  phase  change  induced  by  the  absorbed  10. 6y 
radiation . 

f.2b  Experimental  Results  and  Preliminary  Analysis 

The  data  recordings  presented  in  the  following  figures  were 
taken  with  the  experimental  arrangement  shown  in  Figure  2 
Figure  3  shows  typical  results  for  our  KC1  samples.  These  data 
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were  reproducible  for  different  samples  and  were  insensitive  to 
changes  of  the  coupling  fluid.  The  strip  chart  recording 
shows  the  phase  variation  of  the  surface  wave  vs.  time,  when 
three  successive  10. 6p  laser  pulses  of  1  second  duration 
were  incident  on  the  KC1  sample.  We  note  a  rapid  increase  in 
the  phase  value  of  Ac^,  during  the  time  of  the  laser  pulse,  and 
then  a  decay  to  a  value  A<t>  above  the  initial  level  after  the 
10. 6y  pulse  ends.  This  step-like  behavior  was  consistently 
observed  when  the  sample  was  irradiated  by  successive  10.60 

pulses.  Figures  4  and  5  show  the  variation  of  A<b  and  the 

Ad>  ,  1 

average  slope  ^  with  laser  pulse  energy  at  a  constant 

pulsewidth  of  1  second. 

Figure  6  shows  the  variation  of  the  ratio  A<J>  / A<(>  with 
the  pulse  duration.  The  fact  that  this  ratio  was  found  to 
decrease  to  1  for  longer  pulses  indicates  the  presence  of 
some  "grow-in"  process.  This  could  be  the  result  of  thermal 
diffusion  bringing  the  acoustically  probed  region  to  equilibrium 
with  the  remainder  of  the  sample,  or  the  onset  of  effects  due 
to  the  constrained  state  of  the  irradiated  region.  The  thermal 
spatial  gradiant  caused  by  the  10. 6p  beam  will  induce  stresses 
in  the  sample  which  tend  to  constrain  thermal  expansion  in  the 
heated  region  during  the  laser  pulse.  When  the  pulse  ends, 
diffusion  will  begin  to  reduce  this  gradiant  and  thermal  expansion 
will  proceed.  Since  the  phase  shift  is  sensitive  to  the  acoustic 
path  length,  it  may  be  indicating  an  actual  increase  in  magnitude 
due  to  this  expansion.  These  possibilities  necessitate  the 
derivation  of  a  mori  accurate  relationship  between  the 
temperature  change  and  the  resulting  phase  variation  which 
would  include  the  effects  of  diffusion  and  thermoelastic  stress. 
This  theory  has  been  completed  and  is  being  applied  to  data 
analysis  at  the  present  time. 
f . 2c  Research  Plans 

In  order  to  minimize  the  effects  of  thermal  diffusion,  we 
will  attempt  acoustic  probe  measurements  with  shorter  laser 
pulses,  perhaps  as  short  as  a  few  milliseconds.  In  this  time 
the  diffusion  process  will  tend  to  thermalize  only  material 
within  the  acoustic  energy  volume  which  for  propagation  on  KC1 
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is  ~70y  deep.  This  will  require  laser  powers  as  high  as  50  - 
75  watts,  based  on  sensitivity  results  of  our  current  measure¬ 
ments.  This  source  is  presently  being  constructed. 

We  will  apply  our  technique  to  absorption  measurements  of 
several  samples  grown  at  Hughes -Mai ibu  laboratories.  These  will 
include  RAP-grown  KC1  with  polished  and  HC1  etched  surfaces  and 
also  polycrystalline  samples.  These  experiments  will  indicate 
the  sensitivity  of  these  acoustic  measurements  to  surface 
quality  and  will  supplement  results  reported  by  other  groups  on 
similiir  materials  and  surfaces. 

We  are  also  preparing  an  interdigital  transducer  pair 
which  will  allow  probe  measurements  on  KC1  at  the  first  and 
third  harmonic  RF  frequencies  with  equal  efficiency  during  the 
same  experimental  run.  In  this  way,  the  experimental  parameters 
can  be  kept  identical  for  two  acoustic  probe  depths  differing  by 
a  factor  of  3.  These  experiments  may  indicate  the  presence  of 
differing  bulk  and  surface  properties. 


F.  A.  Kroger,  J.  H.  Margurger  et  al.,  Electronic 
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F19628-72-C-0275,  ARPA  Order  No.  2055. 

ibid.,  Report  No.  3  (December  -  February  1973) 
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ARPA  Order  No.  2055  (29  September  1972). 
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-FIG- 2  General  Experimental  Arrangement  for  Reading 

the  Temperature  Induced  Change  in  Output  Phase,  A<p , 
Brought  About  by  the  Absorption  of  a  C02  Laser  Beam. 
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g.i  Characterization  of  Optical  Performance  of  IB.  Window  Systems 
J.  H.  Marburger ,  Martin  Flannery 

of  th  VT,  llUS  qUarterly  reP°rt'  WC  be84"  a"  extensive  review  of  the  theory 
o  thermal  leasing,  which  summaries  work  performed  on  this  contract,  and 

a  so  tnclndes  results  first  obtamed  by  others,  e.g.  Bendow  and  Cianino  at 
fcRl.  Thls  review,  wh.ch  w,u  be  publjshed  elsewhere  whcn  comple[e_ 

1  S°mewhat  m  viewpoint  from  previous  work,  and  stresses  physical  inter, 
pretation  of  the  formalism, 

l"  The  Diffraction  Integral  for  Beams 

The  starting  point  for  all  diffraction  computations  in  this  work  is  the 

slowly  varying  envelope  equation  for  the  complex  optical  field  envelope  2 

which  is  related  to  the  optical  electric  field  j£  through 

opt.  6 

E  =  Re  Eei(kz  "  a)t) 
opt  ce 


0.1) 


ere  z  is  the  dtrection  of  propagatton  of  the  beam  under  consideration,  and 
-  n»/c,  n  being  the  refractive  index  of  the  propagation  medium  (air)  at  the 
dominant  fluency  „  of  the  beam.  Because  spatial  variations  along  the  beam 
irecfonare  supposed  to  occur  over  distances  large  compared  to  X  =  2„ /k  we 
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Under  these  conditions,  Maxwell's  equations  for  propagation  in  a  linear  isotropic 
homogeneous  medium  imply 
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Here  V  is  the  group  velocity,  and  we  have  ignored  group  velocity  dispersion. 
vt  1S  the  transverse  laplacian 
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With  the  change  of  variable 


t  -  z/v  eqn.  (1.  3)  be 


comes 


where  S  i: 


■  ,  3  b  2-> 

2lk  sT  +  ’te 


now  a  function  of  x,  y,  z  and  r.  S 


(1.4) 


ince  this  equation  does  not  depend 


>n  T,  all  lime  dependence  is  contained  in  the  initial  condition,  whtch  we  write  as 

£0(x,y,t)  =  E (x ,y , o , t ) . 


(1.5) 


The  general  solution  of  (1.  4)  with  initial  condition  (1.  5)  is 


E(r,z) 
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(1.6) 


where  the  t  dependence  is  suppressed,  and  r=(x,y).  If  g  ts  derived 
from  a  fteld  ^  which  passes  through  a  thin  positive  lens  of  focal  length  R  at 
Z  =  o,  then  Eq  may  be  replaced  in  (1.  6)  by 

=  exP  -iJ<r?72R  .  (1>  7-, 


The  incident  beam  ^  is  usually  limited  by  an  aperture  of 


dimension  d,  which  is  a  cc 


characteristic 


.  tonvement  scaling  length  for  the  transverse  variables. 

Ihe  near  field  length  of  the  anerhirp  i „  j  2  -in.-. 

s  cne  aperture,  kd  ,  is  a  convenient  scaling  length  for 

longitudinal  variables.  Introducing 


P  =  r/d, 


C  =  z/kd“,  o  =  R/kd2, 


we  find  that  the  diffraction  integral  is 


(1.8) 
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The  variables  p  and  Q  are  somewhat  more  easily  visualized  than  the 

?  1 

conventional  u,  v  variables  employed  by  Born  and  Wolf,  and  Bendow  et  al. 
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The  diffraction  kernel  K  is  simpler  in  appearance  in  the  variables  u  and  v, 
but  for  complicated  diffraction  fields,  it  is  useful  to  employ  variables  like 
p,  Q  which  are  directly  proportional  to  transverse  and  longitudinal  distances, 
with  constant  scale  factors.  We  use  both  sets  of  variables  to  make  inter¬ 
mediate  formulas  more  concise,  but  all  graphs  and  interpretations  will  involve 
only  p  and  Q  . 


2.  Vector  Corrections 


The  slowly  varying  envelope  equation  (1.3)  is  satisfied  by  field  vectors 
whose  polarization  remains  parallel  to  the  plane  of  incidence.  But  the  trans- 
versality  condition  V  •  ~  0  for  the  exact  optical  field  requires  that  the 

local  plane  of  polarization  be  tangent  to  the  local  wave  front,  or  normal  to  the 
rays.  Thus  z  •  E?  must  be  proportional  to  sina,  where  a  is  the  angle 
between  the  local  ray  direction  and  the  z  axis.  This  is  on  the  average  less 
than  d/R,  the  inverse  of  twice  the  f  number  of  the  window  system.  Thus 
a  =  •  1  for  an  f5  system.  This  angle  must  be  small  for  the  slowly  varying 
envelope  approximation  to  be  accurate,  so  the  error  incurred  in  setting 
z»f!  =  0  is  small.  It  is  nevertheless  useful  sometimes  to  know  the  magnitude 
of  the  longitudinal  field  (e.g.  near  a  focus).  This  is  very  nearly  equal  to 
the  z  component  of  the  projection  of  E?  on  the  plane  normal  to  the  local  ray 
direction  n.  That  is,  the  field 


(I  -  nn)E  , 


is  a  good  approximation  to  the  transverse  optical  field  amplitude  in  the  context 
of  the  slowly  varying  envelope  approximation.  The  vector  correction  -n  (n«E) 
is  ignored  in  what  follows. 
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3.  The  Aberration  Matrix  for  Thin  Windows 

The  initial  field  e!  of  interest  in  window  studies  is  the  field  at  the  exit 
plane  of  the  window.  This  exit  field  may  be  obtained  from  the  field  E?.  incident 
on  the  window  by  solving  the  inhomogeneous  envelope  equation 


2ik  ^  +  V7^p 

i  57  +  vte 


(3.1) 


in  the  window  region  -w  SzSO.  Here  =  n^uu  /c,  and  n^  is  the  mean  scalar 
refractive  index  of  the  window.  This  equation  includes  the  linear,  inhomogeneous, 
possibly  anisotropic  susceptibility  X°f  the  window  system.  If  the  window  thick¬ 
ness  w  is  less  than  the  scale  of  variation  of  X.  and  also  less  than  the  shortest 
diffraction  length  of  any  intensity  or  phase  fluctuation  in  the  incident  beam,  then 
the  approximate  solution  of  (3.1),  ignoring  reflections  at  the  surfaces,  is 

£(r,0)  =  [exp  ^i2TTk1w^ (r) /n|) j  *E(r , -w)  .  (3.2) 


A  (r  )  •  E  (r  ,  -  v) 


The  quantity  X  (~r)  will  be  called  the  aberration  matrix  for  a  thin  window.  If 
any  of  the  assumptions  leading  to  (3.2)  are  violated,  we  have  a  thick  window, 
for  which  the  analysis  is  much  more  complicated  and  will  not  be  treated  in  this 
report.  (The  no- reflection  assumption  is  independent  of  window  thickness  and 
will  be  discussed  later.) 
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The  goal  of  thin  window  theory  is  to  compute  the  2x2  matrix  x(r)  from 
the  known  properties  of  the  incident  beam  and  the  medium,  and  to  determine 
the  influence  of  these  properties  on  the  diffracted  field  S  (r,z)  beyond  the 
window.  Figure  (3.1)  shows  the  geometry  of  our  window  model. 

4.  Aberration  Matrix  for  Cylindrical  Anisotropy  and  Inhomogeneity 

(Thin  Windows) 

An  important  class  of  window  problems  has  aberration  matrices  which 
are  diagonal  in  a  cylindrical  basis: 


X(?) 


A  fr)e  e  +  A  fr)e^e_. 
r v  '  r  r  J  0  0 


(4.1) 
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Fig. 


\ 

I 


Coordinate  system  and  definitions  of  fields.  E  is  the  unfocused  incident  field,  and 
Ej  the  aberrated  field  after  passing  through  the  window.  E  is  the  focused  aberrate 
field  at  z-0  and  E  the  field  at  z>0.  r'  or  r1,  0'  locates  points  in  the  z  =  0  plane  and  r 
°r  1,  0  locates  pgints  in  the  z  plane.  Formulas  are  usually  simpler  in  the  variables 
P-r/d,  Q  -  z/kd  where  d  is  the  window  radius. 


Here  e  and  e  art  radial  and  azimuthal  unit  vectors,  and  A  ,  A  are  the  local 
r  u  r  0  -> 

eigenvalues  of  the  aberration  matrix.  Some  important  properties  of  A  are: 


at  -  J, 


->•  -v-L  -4 

A-A'  =  I 


TrA  =  Ar  +  AQ, 


det  A  =  ArAQ 


(4.2) 


A-A  =  (TrA)X  -  (det  A)I  . 


One  usually  works  in  a  cartesian  basis  if  the  incident  field  is  linearly 
polarized.  In  this  basis. 


t  - 


/A 

A 

I  XX 

xy 

\  A 

A 

yx 

yy 

Ar  cos20  +  Aq  sin20  (Ar-A0)  sin0  cos0 
(Ar”A0)  sin0cos0  Ar sin20  +  AQ cos2© 


(4.3) 


where  r,  0  are  the  polar  coordinates  of  the  point  in  the  window  at  which  A  A 

,  ,  r  6 

are  evaluated. 
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Since  X  is  unitary,  its  eigenvalues  may  be  written,  (in  the  noLation  of 
Bendow  and  Gianino), 


=  exp  ik1$r(rj 


A  =  exp  ik  t  f r) 

0  i  0 


(4.4) 


If  we  introduce  the  mean  and  Qiiference  phases  througli 


^kl  ( 4>r  +  4>0) 


'  V' 


(4.5) 
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we  find  the  following  convenient  form  for  X: 


A  =  e' 


cosA  +  s inA  cos 20 


sinA  sin20 


sinA  sin20 

cosA  -  sinA  cos20 


(4.6) 


id)* 

e  . 


If  the  window  is  isotropic,  A=0,  £  =  f,  and!  =  ?eicP.  (4.7) 

Explicit  formulas  for  A  and  cp  will  be  derived  in  subsequent  sections  for 
specific  window  models.  Cylindrical  anisotropy  generally  requires  that 
all  stresses  and  heat  sources  within  the  window  be  distributed  with  cylindrical 
symmetry  so  that  the  local  principal  polarizations  are  radial  and  azimuthal. 

5.  Diffraction  Field  for  Cylindrical  Anisotropy 

The  aberration  matrix  allows  us  to  find  the  field  E*  in  the  diffraction 
integral  (1.9)  in  terms  of  the  unfocussed  field  £  =  f*  (r,-w)  incident  on  the 
window.  Denoting  the  difference  in  azimuth  between  source  point  and  field 
point  by  =  0'  -  0,  we  may  rewrite  eqn.  (1.9)  as 


where 


E(P,6)  - 


1_  exp  6^]  x 
ir  i  ?  y  2 J 


x  1  p  '  dp  '  exp  (- %iu  p'  2)  exp( i cpjA?*  (  p' )  , 

J  n 


G  =  f  dip  exp  (- ivp  '  cos^) B (p  '  ,0  ' ) 


(5.1) 


JqcosA  -  J2sinAcos20 


^sinA  sin20 


-J^sinA  sin20 


Jq,osA  +  J^sinA cos20  .  (5.2) 


Here  Jq  and  are  Bessel  functions  whose  arguments  are  v  p'  =  po  '/£,  and 
0  is  the  azimuth  of  the  field  point  at  (  p,  Q). 
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Ihe  dependence  of  E  (  p,  C)  on  0  is  explicitly  displayed  in  the  matrix  (5.2), 
To  discuss  this  dependence,  it  is  convenient  to  rewrite  the  diffraction  field  in 
the  concise  form 


E(P,0,C)  =  E0(p,0  -  R(0)-E2(p,?) 


(5.3) 


where 


3(0) 


cos  20  sin  20 

sin  20  -cos  20 


(5.4) 


E  (p,C) 

o 


L-^p2/2cf10.  do' 


2  7T  i  z; 


e  v  cosA  J  E- 

o  1 


(5.5) 


*20.0  -  1ST  e±o2/25/„lp'  dp 


,  -^lUP  2  1$  .  .  ,  t 

e  e  smA  J„E . 

2  i 


(5.6) 


The  matrix  it  is  its  own  inverse,  so  the  intensity  is  proportional  to 


E (p ,0 , C) 


E0|2  *  |t2|2  -  2Re^.R(S).«J 


(5.7) 


This  important  result  shows  that  the  intensity  distribution  beyond  the 
window  is  not  cylindrically  symmetric  for  a  general  incident  polarization. 

The  actual  angular  shape  of  the  intensity  distribution  depends  upon  the  incident 
polarization  state  and  also  upon  the  relative  phases  of  and  £  Notice  that 
for  an  isotropic  window  (no  stress  indi  ced  birefringence)  A  =  0,  £^=0,  and 


E  ( o,  0,  C  ) 


Eo  (  D >  C ) 
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Therefore  the  departure  from  cylindrical  symmetry  in  the  diffraction  field 
IS  caused  by  anisotropy  (e.g.  stress  induced  birefringence). 

6‘  ^M.fracted  Beam  Geometry  for_Specific  Incident  Polarization  States 

Eqn.  (5.  7)  for  the  intensity  beyond  the  window  may  be  evaluated  to 
determine  the  explicit  angular  dependence  for  specific  incident  polarization 
states.  Simple  results  are  obtained  only  if  the  incident  polarization  is 'the 
same  for  all  points  in  the  window  plane,  (e.  g.  linear  polarized  in  same 
irection  at  center  as  in  periphery).  In  this  case  g.  may  be  written 


g. 


i .  f 
l 


ihy 


cosp  e  +  e 


-ihy 


sinn  e 


y. 


(6.1) 


where  E.  is  a  function  of  r'  or  p'  alone,  and  the  quantities  Y, 
to  the  orientation  angle  T  and  ellipticity  e  (ratio  of  minor  to 
Of  the  incident  polarization  ellipse  through 


r\  are  related 
major  axis) 


tan  -  tan2p  cosy 
e  =  tany 

sin  2X  =  sin2n  siny  . 


(6.2) 


Because  the  incident  orientation  is  irr 
The  computation  of  £ 

o  2 

trigonometry,  one  finds 


elevant,  we  shall  take  y  =0  in  the  following, 
is  now  straightforward.  After  some 


Eo-i?(0) 


■k 

fcQE2  cos2X  cos 20 


(6.3) 


where  Eq  and  E2  are  obtained  from  (5.  5)  and  (5.6)  by  replacing  &.  by  E.. 
The  angular  distribution  *  * 


lE0l  +  I  E  2  I  ^  '  2  (Re  FbE2)  cos2x  cos  20  (6.4) 
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n  uiar  plo  olin  tensity  at  fixed  r,z  for  linear  polarization  and 

,,o  2'  ~  *  l.oEzKm  These  curves  resemble  but  are  not  actually 
he  constant  intensity  contours,  which  are  sensitive  to  the  detaile. 
ependence  of  E  ,E  on  radius.  For  the  curves  shown, 

intlnsity°|Eo"|2^'+  |  °f  the  circle  15  the  “"modulated 
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therefore  oscillates  between  the  shapes  plotted  in  Figure  (6.1)  as  the  relative 
phase  of  E^  and  E^  changes.  The  angular  modulation  is  most  severe  for 
linear  incident  polarization  (2y  =  0),  and  vanishes  for  circular  polarization 
(2X  -  ±  tt/2).  Cylindrical  symmetry  is  also  preserved  for  unpolarized 
incident  light,  for  which  the  diffraction  field  is  an  incoherent  superposition 
of  all  possible  polarization  states.  Thus 

2  2  2 

|£|  =  I  E0  |  +  |  E  2 |  (unpolarized,  circularly  polarized)  (6.5) 


+  E. 


2  (Re  E  E_)cos2Q  (linearly 


(6.6) 


polarized) 


7.  Relation  to  Notation  of  Bendow  and  Gianino 


Bendow  and  Gianino* decompose  S  into  cartesian  components,  and 
define  through 

$  =  (?Tri?)  exp(-ip2/2?)  £  =  b2 

to  remove  the  initial  factor  in  (5.1).  The  components  of  6,  defined  through 


t)  =  U.  e  +  IF  e 
lx  2  y 

are  related  to  our  fields  by 

o  2  7 


(0  =  0) 


(7.1) 


b"  1U, 


b-h. 


E  -  E„ 
ox  2x 


E  +  E  0 
oy  2y 


(7.2) 


This  simple  relation  holds  only  because  Bendow  and  Gianino  have  assumed 
that  e^  lies  in  the  plane  containing  the  origin  (r'-z=0),  the  focal  point 
(r  -0,  z-R)  and  the  field  point  (r,  z).  This  means  that  their  formula  for  l5 
is  restricted  to  the  0=0  axis.  A  simple  formula  for  0^0  involving  and 
is  possible  only  lor  unpolarized,  circularly  polarized,  or  linearly  polarized 
light. 

For  linearly  polarized  light  whose  incident  polarization  makes  angle  f 
with  the  x  axis,  the  diffraction  field  envelope  at  (r,  2=0,z)  may  be  written 
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(7.3) 


£  = 


(E  - 
0 

^2 )  cos'! 

(E  + 
o 

E  2  )  s  i  n  V 

where  Eq  and  E.,,  defined  beneath  eqn.  ,6.3),  are  proportionai  to  the  incident 
field  amplitude  along  the  polarization  direction.  In  reference  3  ,  the 

intensity  corresponding  to  (7.  3)  is  written  as 


b|2|g| 


cos2f  |u  .  (i;  ->E  )  | 2 

1  X. 


+  sin2»|U2(B  -E)|2  . 


For  more  complicated  incident  polarization  states,  eqn.  (6.4)  must  be 
employed.  When  the  diffraction  field  is  expressed  in  terms  of  the  fields  E 
an  E2,  the  resulting  formulas  are  usually  simpler  in  appearance,  and  eas°er 
o  understand  from  a  physical  viewpoint  than  corresponding  formulas  based 
upon  Uj  and  U.,.  We  emphasize  that  E.,  vanishes  and  E  becomes  E  when  there 
is  no  birefringence.  Later  we  shall  discover  more  physically  meaningful 
differences  between  E.,  and  E<),  and  use  them  to  develop  our  understanding 
o  tne  influence  of  birefringence  on  the  transmitted  field. 

8.  General  Properties  of  Radial  Dependence  of  Transmitted  Intensity 
A.  Axial  Intensity 

dependence  on  r  or  p  =  r/d  is  contained  in  the  component  amplitudes 

C)  and  F,  crii/on  k.-/c  c  \  _  i  t  r-  /, 


-  '  -  Lrie  component  amplitudes 

th° e  °V  “h  2  °'C  .RiVe"  by<5-5)  and  (5'6)-  In  equations  p  appears  in 

tnG  -/litial  phase  function  \x/Viir-v.  a  ^  ~  ^ <  ... 


.  ,  _  ,  1  equdcions  p  appear 

initial  phase  function,  which  does  not  contribute  to  the  intensity,  and  in 

•he  arguments  of  J.  (p  p7£>  and  p  p7£>.  Since  J  ,0,  =0.  the  field  E  (0,  ,  , 
vanishes,  and  the  on-axis  field  is  ^ 


S(0,  0,  (,)  =  fo(0,  C) 

ThU,S  ~ polarization  as  the  field  incident  on  thf 
^fej.v?n  iil  presence  of  birefr inPenr.  - ~ 
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At  the  geometrical  focal  point,  z-R,  u  0,  the  magnitude  of  E  is 

o 


Eo(0,?) 


rl 

-  p'dp'e1' 
*  Jn 


cosAE^(p' ) 


(8.2) 


Both  cp  and  A  are  functions  of  o'  here,  and  for  cylindrical  symmetry  they 
are  expected  to  be  even  functions.  One  hopes  that  the  phase  of  E.  is 
constant,  or  at  least  a  weak  function  of  p',  so  that  the  phase  variation  of 
the  integrand  is  dominated  by  cp .  If  cp  and  A  are  expanded  in  powers 

of  p'  , 


t  y  ,  i  4 

4>  =  <J>0  +  4>2p  z  +  d4P  +  . . . 


'2  '4 

A  =  Aq  +  A2p  +  A^p  +  ...  , 


(8.3) 


then  cpQ  contributes  a  time  dependent  phase  which  does  not  affect  the  intensity, 
and  cp2  P^  has  the  same  form  as  the  phase  introduced  by  a  thin  lens  of  focal 
length 


-k/2<|>. 


(8.4) 


This  may  be  positive  or  negative  depending  upon  the  sign  of  cp^,  and  may  be 

removed  by  adjusting  the  focal  length  R.  The  remaining  terms  in  cp 

cause  spherical  aberration,  which  reduces  the  focal  intensity  by  an  amount 

depending  upon  the  distribution  of  the  incident  intensity  in  the  window 

aperture.  If  the  aperture  is  filled  uniformly,  then  the  focal  intensity  is 

4 

reduced  appreciably  when  cp^  «tt/2  (the  "Rayleigh  limit").  An  exact  compu¬ 
tation  for  truncated  gaussian  beams  and  gaussian  aberration  is  included  in 
a  later  section. 

The  lowest  order  in  the  expansion  of  A  would  contribute  to  the  inten¬ 
sity  if  it  were  finite,  but  at  the  point  p'  =0  the  window  must  be  isotropic 
because  of  the  cylindrical  symmetry,  and  thus  Ao=0.  Therefore  cosA(p') 
is  unity  at  the  origin,  oscillates  where  A  is  large,  at  some  finite  value  of 
o',  and  then  returns  to  unity  if  the  window  is  isotropic  near  its  pe.1  ’^hery. 
This  has  the  effect  of  removing  annuli  from  the  incident  field  at  the  radius 
where  A  has  its  maximum  value. 
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I  he  overall  result  is  reduction  of  the  focal  intensity  by  a  factor  which  may 

oc-iliate  With  increasing  A.  A  more  detailed  understanding  of  the  influence 

of  A  may  be  gamed  in  the  context  of  particular  models,  developed  in  a  later 
section. 

At  axial  points  other  than  the  focal  point,  we  expect  structure  in  the 
intensity  whose  complexity  increases  as  the  aberration  increases,  and  also 
as  the  transverse  spatial  structure  of  the  incident  field  becomes  more 

abrupt.  For  example,  the  axial  intensity  of  a  uniformly  illuminated 
aperture,  easily  found  from  (8.1)  and  (5.  5),  is 


E  ( 0 , 0  ,  c  ) 


TTU  £ 


i  .21 

‘T  sin  fu 


O  -  O 


2  Sin2?(?  ‘  |)  •  (8.5) 


This  has  an  infinite  number  of  maxima  and  minima  prior  to  the  focal  point 

and  few  beyond  it.  In  contrast,  the  axial  intensity  of  an  untruncated  gauslian 
beam,  of  1/e  intensity  radius  a,  is 


|E (0,0, C') | 


L  /  »  ? 

*  n  (i  +  u  l) 


;T7TF?Tj 


(8.6) 


7  '  Z/ka  ’  R/ka2-  Th'“  “  *  Action  w„h  only  one  maximum 

a  e  beam  wa, st.  The  presence  of  the  abrupt  edge  at  r  d.  z  0,  in  the  field 

t.ie  i  lu.ninulou  aperture  is  clearly  responsible  for  the  violent  nrefocal 

oscillations  in  the  transmitted  intensity  ,8.5).  In  general,  increasing  truncation 
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andjninima  of  increasing  amplitude  prior  to  the  focal  region. 


■B.  Off-Axis  Intensity 


Equations  (5.5)  and  (5.6)  for  the  fields  Eq  and  E2  are  simplest  in  the 
geometrical  focal  plane  z  =  R  or  Q  =  $ : 


eip2/2  f1  id) 

2 tt i^1  J  P'dP’e  ?cosA  J  (pp' /  ^ )  Ei  (p- )  , 


(8.7) 


E  2  (p  J?  ) 


eip2/2  f1  id, 

“Hi?  J  P'dP'e  9  sinA  J2(ppVp)Ei(p') 


(8.8) 


Because  J2  vanishes  at  zerc  argument,  E0  contributes  only  to  the  off-axis 
intensity  In  particular,  the  intensity  contribution  of  E2  at  radius  p  is  most 
strongly  influenced  by  the  value  of  FsE.  sinA  expi  cp  at  p'  w  3. 1  f/p  =  3.  1 
R/kdr,  where  has  its  first  maximum  (J2(3.1)  «  .486).  Now  A  vanishes  at 

p1  =0,  30  |f|  possesses  a  maximum  at  some  finite  value  p'=p'  .  Therefore, 

!E2  j  is  exPected  to  be  largest  near  pm  =  3.1?/p^.  In  many  ca^es  A  is 
initially  zero  when  the  beam  first  strikes  the  window,  and  then  increases  to 
a  thermal  equilibrium  value.  As  A  increases  we  expect  to  see  a  ring  form 
in  the  focal  plane  near  p^.  The  intensity  of  the  ring  may  oscillate  if  A(p'  ) 
increases  beyond  n/2.  m 

The  component  Eq  tends  to  have  its  maximum  at  P  =  0  when  cp  and  A 
vanish,  because  Jq  has  its  maximum  value  at  zero  argument.  The  effect 
of  truncation  and  aberration  is  as  usual  to  spread  the  energy  over  a  broader 
central  maximum,  and  to  create  subsidiary  rings  in  the  focal  plane. 

In  conclusion,  we  may  view  the  effects  of  cp  and  A  on  the  focal  plane 
intensity  pattern  as  follows:  The  initial  focal  point  is  near  z  =  R,  but  drifts 
to  larger  or  smaller  values  of  z,  depending  upon  the  sign  of  d 2  cp  /dr'2,  as 
the  mean  aberration  cp  increases.  At  the  same  time  the  axial  intensity  peak 
becomes  broader  and  weaker  and  off  axis  maxima  begin  to  appear.  The 
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111  '"’6Uy  "  dePart“«  cylindrical  symmetry  of  this  off-axis  structure 

tncreases  with  the  induced  birefringence.  The  intensnty  pattern  near  the  axis 
where  EQ  donninates,  is  neariy  cylindrically  symmetric,  while  the  ring  asso-  ' 
ciated  with  E2  has  angular  maxima  or  minima  along  directions  parallel  to 

e  maJOr  iX1S  °f  the  incident  fcltnrixat.on  ellipse.  Oscillations  in  the  ring 
intensity  and  in  the  axial  intensity  tend  to  be  out  of  phase  because  E  is 


proportional  to  the  integral  of  cosi,  and  E  to  the  integral  of  sinA. 


o 


9'  ^^-r--ation  Functions  for  a  Simple  Window  Model 


To  increase  our  understanding  of  the  aberration  functions  cp  and  A  , 


let  us  commute  them  for  the  model  analysed  extensively  by  Sparks5.  Jasperse 

and  Gianino,  and  Bendow  et  al.  In  *-his  model  •  j 

nis  model,  the  window  is  assumed  to  be 


eated  by  a  cylindrically  symmetric  gaussian  beam,  and  the  local  temperature 

“ hSC  a7  °nly  by  fte  heatin6  rate'  affected  by  thermal  diffusion. 

e  window,  which  is  isotropic,  expands  and  deyelops  internal  strains  as 
he  cooler  periphery  restricts  the  expansion  of  the  warmer  interior.  It  is 
ree  y  supported  so  the  stress  vanishes  at  the  window  boundary.  The  two 


--  vtu,  x  )  .  -L  lie  L WO 

berration  functions  ^and  for  this  case  have  been  derived  by  the  authors 
men  tone  a  ove,  and  presented  in  convenient  form  by  Bendow  and  Gianino1; 


$ 


=  C"f(p) 


n 


C?g(p) 


+  C 


n 


n 


P  ,9 


0.1] 


Here  f  and  g  are  functions  which  depend  only  upon  the  incident  beam 
pa  rameters : 


f  Cp) 
e  Cp) 


_  ?  ? 
exp  -2a  “p 


1  -  exp  (-2a 


V)J/a2p2 


(9.2) 

(9.3) 


where  2a2n^  -  / 2  .  ,  .,  ,  , 

/a  ,  with  a  the  1/e  intensity  radius  of  the  incident  gaussian. 
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/  ,s  a  truncation  parameter.  The  constant  phase  correction 

n  u  kn  - - .*  ri 


=c«  has  no  infiuenoe  on  the  intensity  and  is  usual, /supp  "“7 The 

remaining  constants  are  functions  of  the  paramete-s  of  the  medium,  and  the 

time  integral  Hit,  of  the  axial  intensity  multiplied  by  the  window  thickness 
at  I  an  absorption  constant: 
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H  (t )  =  CBw/C’ ) /tI0Ct  ’  )dt* 


(9.4) 


Here  Sn/dT  is  the  rate  of  change  of  index  with  temperature  at  zero  stress, 
a  is  the  coefficient  of  thermal  expansion,  v  Poisson's  ratio,  and  p  the 
elasto-optic  tensor  elements,  g  is  the  absorption  coefficient  and  C1  the 
specific  heat. 

With  this  notation,  our  aberration  functions  are 


4>  =  +  ) 

=  ^(Cf  +  C®)f(p) 

A  =  k  x  C  2  [g  ( P )  -  2  f  C  P )] 

Thus  in  this  formulation  each  aberration  function  is  proportional  to  one 
lumped  material  parameter,  and  a  function  of  c  depending  only  upon  the 
incident  field.  If  the  incident  intensity  is  constant  in  time,  then  H(t)  is 
proportional  to  time,  and  cp  and  A  are  proportional  to  the  dimensionless 
time  t  : 


(9.5) 

(9.6) 


<J>  =  i  f  (  p  ) 


(9.7) 
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w  h  e  r  e 


A  =  or  [g(p)  -  2f(p)] 
2 


(9.8) 
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It  is  important  to  notice  that  in  this  formulation,  the  introduction  of 
birefringence,  /  p^,  requires  the  introduction  of  only  one  additional 

material  parameter  g  into  the  equations  in  the  context  of  the  model  most 
frequently  used  to  estimate  the  influence  of  birefringence.  The  value  of  a 
may  be^inferred  easily  from  the  ratios  C  °/ C°  ,  C^0/  C®  tabulated  by  Bendow  and 
Gianino,  and  is  shown  in  Table  \9.1)  tor  the  materials  they  considered,. 

Now  let  us  examine  the  functional  dependence  of  A  and  cp  upon  p. 
Figure  (9.  1)  shows  f  and  |(g-2f)  plotted  on  the  same  scale.  The  function 
A  /a  T  reaches  its  maxir  am  value  of  .  298  at  r/a  =  1.  34  where  0  /T  has 
declined  to  .167.  This  means  that  the  birefringence  is  largest  in  the  low 
intensity  wings . 

With  these  expressions  for  the  aberration  functions  cp  and  A,  we 
shall  analyse  the  following  problem:  given  a  particular  incident  beam 
truncation  parameter,  for  what  value  of  the  material  parameter  a  does 
birefringence  begin  to  play  an  important  role  in  reducing  the  focal  intensity? 
To  this  end,  we  first  determine  the  focal  intensity  without  birefringence. 

10-  Focal  Intensity  with  No  Birefringence  for  Simple  Window  Model 

Taking  A=0,  we  find  the  focal  field 
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(10.1) 
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Table  9.1  Dimensionless  Times  and  Birefringence  Parameters 
(Based  on  values  tabulated  in  Ref.  1) 


Material 


r/wj  I  (t1 


KBr 

-  .598 

-  0.  14l. 

KCA 

-1.  344 

-  2.  51 

NaC£ 

-0.  202 

-  6.  30 

KI 

-0.296 

-  7.91 

Csl 

0? 

-  63.6 

CsBr 

0? 

r- 

< 

o 

00 

1 

CdTe 

.0317 

-181. 4 

KRS-5 

0? 

588. 

ZnSe 

.0746 

1424. 

GaAs 

-  .00886 

2990. 

InSb 

0? 

9790. 

Ge 

.0139 

8570. 

TI  Glass^l  173 

0? 

2890. 

TI  Glass#  20 

0? 

3780. 

Irtran-4  (Zr.Se) 

.  0829 

4310.. 

Irtran-6  (CdTe) 

.0423 

19500. 

w  =  window  thickness  (cm) 

2 

I  =  incident  axial  intensity  (Watt/cm  ) 
o 

t  =  time  since  turn-on  (sec.) 

<j  and  t  are  dimensionless 


(10.2) 


where  the  incident  field  is 

E.(o')  =  E.  exp  -  a2p'2 

7 

As  we  have  shown  elsewhere,  this  integral  can  be  evaluated  exactly  in  terms 
of  Fresnel  functions,  which  can  be  easily  visualized  with  the  aid  of  Cornu's 
spira'.  The  result  for  the  focal  intensity  is 
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(10.3) 


2 

Here  h  =  exp  -2a  ,  and  S(u),  C(u)  are  the  Fresnel  sine  and  cosine  integrals 
which  give  the  y  and  x  coordinates,  respectively,  of  a  point  on  the  Cornu 
spiral.  The  argument  u  measures  the  arc  length  along  the  spiral  from  the 
origin  to  the  point  [C(u),  S(u)]. 
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In  the  absence  of  thermal  distortion  (  t  =  0)  the  focal  intensity  becomes 

(10.4) 


PQkV(l  -  J)2 


TT  R‘ 


where  Pq  is  the  total  power  in  the  incident  gaussian  beam,  and  a  is  the  1/e 
intensity  radius.  This  function  has  a  maximum,  if  regarded  as  a  function  of 
a  for  fixed  d,  when  h=  0.0805,  or  a  =  0.63d.  That  is,  this  choice  of  truncation 
maximizes  the  focal  intensity  for  fixed  beam  power  and  window  size,  and 
provides  a  useful  reference  point  for  the  truncation  parameter. 

For  the  optimum  value  h  =  0.0805,  evaluation  of  eqn.  (10.3)  shows 
that  the  focal  intensity  drops  to  half  its  aberrationless  value  when  t=Ti  =0.965  n. 
71ns  characteristic  degradation  time  is  very  close  to  that  estimated  by2Jasperse 
and  Giamno,  who  analysed  this  model  using  geometrical  optics.  We  shall  set 
1I  tt  for  convenience  in  subsequent  sections. 
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he  remainder  of  this  analysis  will  be  included  in  the  next  qua 

T"  1 _ ......  * 
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report.  Ihe  new  results  reported  here  may  be  summarised  as  follows: 

I  he  analysis  of  thermal  lensing  based  upon  the  elementary  fields  g  and  g 
leads  to  simple  formulas  for  the  azimuthal  dependence  of  the  transited  2 
field  intensity  and  polarization  state  for  an  arbitrary  incident  polarization 
state.  It  also  dearly  separates  the  effect  of  birefringence  from  the  "scalar" 
refocusing  and  aberration  phenomena.  Furthermore,  it  allows  the  influence 
of  birefringence  to  be  characterized  by  a  single  material  parameter  o  rather 
than  the  two  parameters  C°/Cf  and  cf/C,*  necessary  in  the  formulation 
of  Bendow  and  Gianino.  This  allows  vast  simplification  both  in  the  develop, 
ment  of  the  theory,  and  in  the  presentation  of  data,  ft  also  provides  an 

intuitive  basis  for  the  much  more  complicated  non-cylindrically  symmetric 
theory  to  be  presented  in  future  reports. 
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DISCUSSION 


The  Wet  absorption  at  10.6pm  found  „lth  melt  gr„„„  CaAs  lemalns 
at  8  ‘  10  cm-  .  It  has  now  been  established  that  the  weak  temperature 
dependence  of  this  absorption  fits  the  predictions  of  theorv  f„r  multl 
phonon  absorption,  supporting  the  view  that  the  mentioned  value  represents 
the  theoretical  limit.  However,  this  conclusion  Is  not  completely  final 
Heak  temperature  dependencies  may  also  result  from  other  mechanisms.  An 
Increase  of  the  IR  absorption  spectrum  after  heating  at  400"C  with 
maintenance  of  the  wavelength  dependence  seems  to  indicate  that  defects 
are  still  Involved.  Attempts  to  grow  GaAs  by  a  completely  different 
method  -  chemical  vapor  transport  -  have  been  partly  successful.  Layers 
up  to  2mm  In  thickness  were  obtained.  Unfortunately ,  the  layers  contain 
Inclusions,  and  the  relatively  high  10. bum  absorption  measured  with  these 
layers  Is  probably  due  to  these.  In  the  next  contract  period  we  shall 
concentrate  on  growing  from  Gad);  these  will  probably  be  of  better  quality 

and  will  enable  us  to  finally  decide  whether  GaAs  Is  a  serious  contender  as 
a  window  material. 

An  electrochemical  method  was  developed  to  measure  oxygen  In  Gad). 

By  dissolving  GaAs  In  this  liquid,  the  oxygen  content  of  the  GaAs  can 

be  determined.  Oxygen  at  a  level  of  10‘'<,-3  was  found  In  a  le-doped 

crystal.  No  0  was  found  In  an  Sl-doped  crystal,  grown  In  oxygen  rich 

environment.  It  Is  believed  that  this  Is  due  to  the  stability  of  S10  . 

The  oxygen  present  In  a  concentration  <  2  x  that  of  SI  Is  not  detected  by  our 
method. 

Structural  Investigations  on  Silicon-doped  GaAs  Indicate  that  hardening 

occurs  In  the  presence  of  tetragonal  defects.  Deformation  studies  show  that 

slip  In  GaAs  occurs  on  only  „ne  system:  (111)  <n0>.  Diffraction  patterns 

show  diffuse  streaks  which  are  attributed  to  diffuse  scattering  of  electrons 
at  acoustic  phonons. 

In  repeating  dielectric  constant  measurements.  It  was  found  that  a 

considera,  le  Inaccuracy  arises  from  errors  in  the  measurement  of  the  sample 

thickness  and  area,  and  In  lack  of  flatness.  The  best  value  found  so  far 

is  k-  13.6,  larger  than  the  value  first  reported.  Further  refinements  are 
being  made. 

Crystals  of  KBr  were  grown  which  showed  an  absorntion  at  10. 6pm  of 
1.1  x  10  3cm  far  above  what  theory  predicts.  KBr03  has  been  suggested  as 
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the  possible  cause  of  this  absorption. 

Attempts  to  prepare  dense  polycrystalline  ingots  of  CdTe  by  hot- 
pressing  in  the  presence  of  an  excess  of  cadmium  still  are  not  completely 
successful . 

We  are  presently  remeasuring  the  high-temperature  Hall  effect  on 
CdTe,  and  will  compare  the  result  with  those  of  measurements  made  on  cooled 
cyrstals  in  order  to  see  what  changes,  if  any,  occur  during  cooling. 
Comparison  of  our  high-temperature  data  on  indium-doped  CdTe  and  published 
data  on  cooled  crystals  indicate  that  native  donors  (probably  interstitial 
Cd)  precipitate. 

Multiphonon  absorption  may  be  due  to  anharmonic  forces  and/or  to  a 
non-linear  dependence  of  ionic  dipoles  or  in  displacement.  An  experimental 
method  to  determine  the  ratio  in  which  the  two  mechanisms  contribute  is 
proposed.  The  theory  of  thermal  lensing  is  being  reviewed  and  is  being 
reformulated  based  on  considering  the  effect  from  a  different  viewpoint. 
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